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Abstract 



Transport properties of magnetic multilayers in current perpendicular to the plane of the 
layers geometry are defined by diffusive scattering in the bulk of the layers, and diffu- 
sive and ballistic scattering across the interfaces. Due to the short screening length in 
metals, layer by layer treatment of the multilayers is possible, when the macroscopic 
transport equation (Boltzmann equation or diffusion equation) is solved in each layer, 
while the details of the interface scattering are taken into account via the boundary con- 
ditions. Embedding the ballistic and diffusive interface scattering in the framework of 
the diffusive scattering in the bulk of the layers is the unifying idea behind the problems 
addressed in this work: the problem of finding the interface resistance in the multilay- 
ered structures and the problem of current-induced magnetization switching. 

In order to find the interface resistance, the method of solving the semiclassical 
linearized Boltzmann equation in CPP geometry is developed, allowing one to obtain 
the equations for the chemical potential everywhere in the multilayers in the presence 
of specular and diffuse scattering at the interfaces, and diffuse scattering in the bulk 
of the layers. The variation of the chemical potential within a mean-free path of the 
interfaces leads to a breakdown of the resistors-in-series model which is currently used 
to analyze experimental data. While the resistance of the whole system is found by 
adding resistances due to the bulk of the layers and resistances due to the interfaces, 
the interface resistances are not independent of the properties of the bulk of the layers, 
particularly, of the ratio of the layer thickness to the mean-free path in this layer. 

A mechanism of the magnetization switching that is driven by spin-polarized cur- 
rent is studied in noncoUinear magnetic multilayers. Even though the transfer of the 
spin angular momentum between current carriers and local moments occurs near the 
interface of the magnetic layers, in order to determine the magnitude of the transfer one 
should calculate the spin transport properties far beyond the interfacial regions. Due 
to the presence of long longitudinal spin-diffusion lengths, the longitudinal and trans- 
verse components of the spin accumulations become intertwined from one layer to the 
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next, leading to a significant amplification of the spin-torque with respect to the treat- 
ments that concentrate on the transport at the interface only, i.e., those that only consider 
the contribution to the torque from the bare current and neglect that arising from spin 
accumulation. 
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Chapter 1 
Introduction 



1.1 Magnetoresistive effects 

During the last decade, attention has been focused on a question of magnetically con- 
trolled electrical transport. Materials can change their resistance in response to a mag- 
netic field. This phenomena is called magnetoresistance (MR). All metals have an inher- 
ent, though small, MR owing to the Lorentz force that a magnetic field exerts on moving 
electrons. Metallic alloys containing magnetic atoms can have an enhanced MR. For ex- 
ample, an anisotropic magnetoresistance (AMR) measures the change in resistance as 
the direction of the magnetization changes relative to the direction of the electric cur- 
rent. (Fig. 11.11) If 9 is the angle between the magnetization and the current direction, 
(Fig. II. lb )), the resistance is the following function of 9: 

R = Ro + ARamrcos^9, (1.1) 

where i?o is the resistance at = n/2, and ARambJ Ro '^^ the AMR ratio. The resistivity 
is typically smaller if the current direction is perpendicular to the direction of magne- 
tization than in the condition that those are parallel due to the scattering anisotropy of 
electrons. The AMR ratio is fairly small: a few percent for Nio.8Feo.2 alloy (permal- 
loy) at room temperature, and somewhat larger at lower temperatures. Nevertheless, the 
phenomenon of anisotropic magnetoresistance has a significant importance for techni- 
cal applications, such as magnetic sensors. As illustrated in Fig. II. lb ), if a magnet is 
attached on a rotating disk, an MR sensor can detect the number of rotations or the speed 
of motion from the resistance change of the MR sensor. It has also been attempted to 
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Figure 1.1: Anisotropic magnetoresistance (AMR): a) Resistance as a function of ap- 
plied field, b) Measurement. V, 1, M, H and 9 are voltage, current, magnetization, 
applied field, and the angle between the current and the magnetization, c) Application 
of AMR sensor. Adapted from Ref. yj 



apply an AMR sensor for a magnetic recording technology, but for ultra-high density 
recording, a very high sensitivity of reading head is necessary, and thus the high MR 
ratio is required, which can't be provided by the AMR effect. 

Magnetoresistive effects are more pronounced in magnetic layered structures, or 
magnetic superlattices. They are formed artificially by alternately depositing on a sub- 
strate several atomic layers of one element, say, iron, followed by layers of another el- 
ement, such as chromium. A term Giant Magnetoresistance (GMR) is used to describe 
the behavior of materials consisting of alternating layers of ferromagnetic and nonmag- 
netic metals deposited on an insulating substrate. The GMR effect has been observed in 
1988 in the resistivity measurements on Fe/Cr multilayers, [2| as shown in Fig. ll.2[ At 
4.2 K the resistivity of the Fe/Cr multilayer was decreased by almost 50% by applying 
an external field. At 300 K, the decrease of resistivity reaches 17%, which is signifi- 
cantly larger than MR changes caused by the AMR effect. If 9 is the angle between the 
magnetization directions of the neighboring ferromagnetic layers, the resistance of the 
system takes the following form: 

R = Rp + ^^^^l-cos9), (1.2) 

where Rp is the resistance of the system when the magnetic moments in the alternating 
layers are aligned in the same direction, or parallel (9 = 0), and ARgmr = Rap — Rp, 
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Figure 1.2: Magnetoresistance of three Fe/Cr superlattices at 4.2 K. The current and the 
applied field are along the same [110] axis in the plane of the layers. Adapted from 
Ref. 

where Rap is the resistance of the system when the magnetic moments of the neighbor- 
ing layers are oppositely aligned, or antiparallel (9 = n). Resistance is the greatest for 
the antiparallel (AP) configuration of the magnetizations, and smallest for the parallel 
(P) one (Fig.Q. The GMR ratio is defined as {Rap - Rp)/Rap * 100% or, alterna- 
tively, as {Rap — Rp)/Rp * 100%. GMR effects has been obtained in two geometries 
(Fig. 11.41) . In the first one the current is applied in the plane of the layer (Current In 
Plane, or CIP geometry), while in the second one the current flows perpendicular to the 
plane of the layers (Current Perpendicular to the Plane, or CPP geometry). The CPP- 
MR is larger than the CIP-MR, and it exists at much larger thicknesses of the samples. 

The mechanism of GMR is attributed to the change of the magnetic structure in- 
duced by an external field. The role of external magnetic field is to change the internal 
magnetic configuration; in cases when it is not possible, the GMR does not appear. It 
is necessary to separate magnetic regions from one another so as to be able to reorient 
their magnetization, otherwise layers are too strongly coupled and ordinary fields can 
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Figure 1.3: Ferromagnetic and antiferromagnetic configuration of a multilayered mag- 
netic metallic structure 
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Figure 1.4: Current In Plane of the layers and Current Perpendicular to the Plane of the 
layers geometries. Adapted from Ref |^ 
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not rotate magnetic moments in one layer relative to moments in the neighboring layer. 
The multilayers are designed so that at zero external field the magnetizations of the 
alternating magnetic layers are aligned opposite to each other. The external magnetic 
field rotates the magnetization, so that all ferromagnetic layers become aligned in the 
same direction, and resistance changes. There are several ways to provide the original 
antiparallel alignment of the neighboring FM layers. The first one involves choosing the 
metallic interlayer thickness (r^ 1 nm) such that the Ruderman-Kittel-Kasuya-Yoshida 
(RKKY) coupling between localized moments via the conduction electrons forces the 
magnetic layers to be antiferromagnetically coupled. The other is to make the succes- 
sive magnetic layers with different coercivities. Smaller field may switch a layer with the 
smaller coercivity, while the layer with a bigger one remains in the previous direction; 
it will rotate at larger field, and there is a range of external fields when the alignment 
is antiparallel. Finally, a phenomena known as "exchange bias" may be used, when an 
AF metal in a contact with FM metal pins the magnetization of the later in the direction 
opposite to the magnetization of the last layer of AF so as to produce zero net magneti- 
zation. The second FM layer, separated from the pinned one by the nonmagnetic layer 
remains free to rotate its magnetization under the influence of the external field. 

The primary source of GMR is spin-dependent scattering of conducting electrons, or, 
more precisely, change in the scattering rate as the magnetic configuration changes by an 
external magnetic field. In the antiparallel magnetic structure, conduction electrons are 
much more scattered that in the parallel magnetic structure. Spin-dependent scattering 
can make a large contribution to the resistivity. In the following section, I will discuss a 
phenomenon of spin-dependent scattering and its effect on the electron transport. 

Similarly to AMR sensors, GMR sensors can also be used to detect rotations and to 
measure the speed of motion. They have an advantage of the full angular dependence. 
Whereas for an AMR-type sensor opposite field directions produce the same signal, for 
a OMR-type sensor parallel and antiparallel magnetization alignments yield different 
resistivities. High sensitivity of GMR sensors make them attractive as the sensors for 
magnetic fields, particularly, as read-out heads in hard disk drives in computers. The 
fact that GMR is mainly an interface effect (see below) allows the sensor to be made 
thinner, which leads to the improvement of the spatial resolution in read-out. GMR- 
based hard drives are already being used in computer industry. The magnetic random 
access memories (MRAMs) based on GMR effect have also been realized, but they are 
not favorable because of the relatively small resistance of the GMR memory element 
compared with the current leads connecting it with the processing unit. 

Besides metallic multilayer systems, similar phenomena are found in granular sys- 
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terns where small ferromagnetic clusters are dispersed in non-magnetic matrices. Ex- 
tremely large MR effect has been found in manganese perovskite oxides and is called 
the colossal MR (CMR) effect. Another class of structures where high MR ratios has 
been obtained is the tunneling junctions, where the layers of ferromagnetic metal are 
separated by a thin insulating barrier; the MR effect in these structures is called tunnel 
magnetoresistance (TMR). 

1.2 Spin-dependent transport in metallic magnetic mul- 
tilayers and the GMR 

The transport in the GMR devices studied to date is dijfusive. In a macroscopic sample, 
electron undergoes scattering off a large number of impurities. One is in the regime of 
diffusive transport when the concentration of impurities requires the averaging of the 
scattering potential. The averaging leads to the loss of the memory of the electron mo- 
mentum direction and, hence, to a resistance. The characteristic lengthscale over which 
electrons retain a memory of their momentum is called the mean free path, A^/p- Due 
to the large transverse size of the layered systems, the macroscopic transport equations, 
such as the Boltzmann equation or diffusion equation, see Sec. El can be applied in 
order to describe electron transport in multilayers, even if the thickness of the layers 
is of the order or smaller than the mean free path. In the multilayers, the presence of 
interfaces leads to additional scattering. At the interface between layers of different 
metals, electrons experience specular scattering due to the band mismatch. The direc- 
tion of the electron momentum changes, but it does so in a predictable and reproducible 
fashion. Specular scattering at the interface between two metals is treated in detail in 
Appendix lA. 3 1 Roughness of the interfaces leads to an additional, dijfusive scattering, 
when the incoming electron unpredictably changes the direction of its momentum, and 
the information about the electron momentum is lost. Due to the short screening length 
in metals (typically about 1-3 A, which is at least an order of magnitude smaller than the 
layer thicknesses in the structures currently being investigated), a layer by layer treat- 
ment of the multilayers is possible, where the macroscopic transport equation is solved 
in each layer, while the details of the interface scattering are taken into account via the 
boundary conditions. 

In ferromagnetic metals, the transport properties of electrons are spin-dependent. 
This dependency arises from the unbalance of the spin populations at the Fermi level due 
to the splitting between the up and down spin states (exchange splitting). As illustrated 
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Figure 1.5: Density of states N{E) for up and down s and p electrons for cobalt, ep is 
the Fermi energy. 

in Fig. ll.5l for cobalt, the majority (spin parallel to the magnetization, or spin-up) d-states 
are all filled, and the d-electron states at the Fermi level contain entirely minority (spin 
anti-parallel to the magnetization, or spin-down) electrons. Although there are also s and 
p electrons at the Fermi level, a significant number of the carriers are the more highly 
polarized d electrons, which produces a current which is partially spin-polarized. The 
imbalance of spin-up and spin-down electrons also results in the different probabilities 
of the s^d transition for majority and minority electrons and consequently in different 
resistivities. 

Spin-dependent scattering occurs both in the bulk of the layers and at the interfaces. 
In a multilayer, electrons move in the potential which reflects the band mismatch at 
the interfaces between the magnetic and nonmagnetic layers. The exchange splitting of 
the up- and down d-bands in the ferromagnetic layers results in different heights of the 
steps seen by up- and down- conduction electrons (Fig. 11.6b )). In the P configuration, 
the height of the steps is the same in all layers but different for majority and minority 
spins. In the AP state, small and large steps alternate for each spin direction. Another 
contribution to the scattering is due to the presence of impurities in the layers, and due to 
the roughness at the interfaces. This scattering is also spin-dependent in ferromagnetic 
metals, and it is stronger for minority electrons (Fig. 11.6b )). resulting in the different 
resistance in P and AP states and the GMR effect. 

A remarkable property of the electron transport in magnetic multilayers is that at 
low temperatures most of the scattering encountered by the electrons does not flip their 
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Antiparallel Parallel 



Figure 1.6: a) Spin dependent potentials in a magnetic multilayer for the AP and P con- 
figurations for majority and minority electrons. Adapted from Ref. [|3l. b) Illustration 
of spin-dependent electron scattering in GMR multilayers. Adapted from Ref. H. 
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spin, as it costs energy. Even though the carrier may undergo many scattering events, 
the orientation of its spin can be very long-lived. The characteristic time when the elec- 
tron remembers its spin is the spin-flip time Tsf, and the characteristic length scale is 

the spin-diffusion length Xg^i = \J Drgf, where D is the diffusion constant. The spin- 
diffusion length is usually of the order of ten times larger then the electron mean-free 
path. The slow spin relaxation leads to a formation of a steady-state non-equilibrium 
build-up of magnetic moments, or spin-accumulation, extending over a length Xsdi from 
the interface between the layers with antiparallel magnetization directions. This magne- 
tization acts as a bottleneck for spin transport across the interface, which in turn hinders 
the flow of charge and results in the higher resistance in the antiparallel configuration 
compared with the parallel configuration. [5 1 

The large spin-diffusion length in ferromagnetic metals means that the current can be 
considered to be carried by two independent channels of carriers, one for spin-up elec- 
trons with a resistivity p|, and the other made of spin-down electrons with a resistivity 
Pj; this is commonly referred to as the two current model. This model allows a simple 
illustration of how the GMR effect works in CPP structures. If the thickness of the lay- 
ers is larger then the electron mean free path in the layers, a resistor network analogy 
is appropriate (Fig. 11.7b )). when the resistivities of each layer for each spin direction 
are added in series, while those for two channels are added in parallel. Since the resis- 
tances of the majority {Rm) and minority {Rm) electrons are different, and Rm < Rm, 
there is a "short-circuit" effect in the parallel configuration, and resistance is smaller in 
P state than in the AP state, and the GMR effect exists. In the CIP geometry, a simple 
resistor network analogy yields the same resistances in both P and AP configurations 
(Fig. 11.7b )). and the absence of the GMR effect. For the CIP-MR to exist, the mean- 
free-path has to be larger than the thickness of a layer. In this case, electrons, which 
travel parallel to the interfaces, sample the scattering from different layers, and the cur- 
rent is sensitive to the change of the magnetic configuration of the system. The electron 
mean-free path is the characteristic length scale for spin-dependent transport in the CIP 
geometry. In the CPP geometry, where the electric current samples scattering in all lay- 
ers, the GMR does not depend on Am/p- The characteristic length scale for the transport 
in the CPP geometry is the spin-diffusion length Xsdi- Spin-flips limit the distance over 
which the two spin currents are independent. If the size of a sample is larger than Xsdi, 
spin-up and spin-down currents are mixed, and the CPP-MR is decreased. In this work, 
I will only consider the current perpendicular to the plane of the layers geometry. 

A CPP-MR model that provides a significant insight into the problem was developed 
by Valet and Fert. Q Within this model, the Boltzmann equation, with an additional 
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Figure 1.7: Two-current model: resistor network analogy for the CIP and CPP resis- 
tances of the ferromagnetic (P) and antiferromagnetic (AP) configurations of the mag- 
netic moments in a multilayer. Rm (red), Rm (green), and Rnm (blue) stand for the 
resistances of the magnetic layers for spin parallel and antiparallel to the local magneti- 
zation, and of the non-magnetic layers. Adapted from Ref. [3J 
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term to account for spin accumulation, is solved at zero temperature. In the case when 
the mean free path is much smaller then the spin diffusion length, the macroscopic trans- 
port equations are obtained that relate the current and the electrochemical potential, so 
that the resistance of the system can be found. The Valet-Fert model correctly describes 
the main features of magnetoresistive multilayers, i.e., increase in MR with decreas- 
ing temperature, increase with decreasing film thickness and increase with increasing 
number of layers. 

1.3 Current-induced switching 

It has been seen from the discussion of the GMR effect that the relative orientation of the 
magnetic moments of the layers affects the electric current, causing different resistances 
for different magnetic configurations. The reverse effect, that a spin-polarized current 
could affect the magnetic moment of a layer, has also been predicted, E [HI |9l [TOl [TTl 
and experimentally demonstrated. iTlllTBlfm In the perpendicular transport geometry, 
the spin-polarized currents may transfer angular momentum between the layers, result- 
ing in the current-driven excitations in magnetic multilayers: either reversal of layer 
magnetization, or generation of spin- waves. [ 15 1 

Spin-polarized current exerts a torque on a ferromagnetic layer if that layer's mo- 
ment is not coUinear with the direction of current polarization. The reversal of the 
magnetization is due to the interaction between the magnetization and the spin accumu- 
lation in a direction perpendicular to the magnetization. Fig. II. 81 shows schematically 
the five-layered structure used to measure current-induced switching, and the directions 
of the torque acting on the magnetic moment of the thin ferromagnetic layer due to the 
spin transfer by the current, polarized in the direction of the magnetization of the thick 
ferromagnetic layer. lITlll 

The great interest in the phenomenon of the current-driven excitations in magnetic 
multilayers lies both in trying to understand the underlying physics and in its poten- 
tial for device use: magnetization reversal for magnetic media and magnetic memories, 
and spin-wave generation for production of high frequency radiation. In present mag- 
netic devices the moments are reversed via externally generated magnetic fields. The 
reading/writing processes would be simplified by applying a polarized current through 
the magnetic layer itself. Practically, though, the magnetization reversal by spin transfer 
requires high current densities, around 10^ A/cm^, in order to overcome magnetic damp- 
ing. The current density should be reduced by approximately an order of magnitude so 
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Figure 1.8: a) Schematic representation of the five-layered structure used to measure 
current-induced switching. The system consists of the thin Co layer of the variable 
thickness (#1), separated by the 4 nm Cu layer from the the thick (100 nm) Co layer 
(#2). Another Cu layer contacts the thin Co layer via an opening 5 to 10 nm in diameter 
in an insulating silicon nitride membrane. b)-e) Directions of torque (red arrows) on the 
magnetic moments in the thin FM layer due to the spin-polarized current. Adapted from 
Ref. yjl 
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that this effect could be considered for applications. 1 16| 

1.4 Problems to be solved 

Transport properties of magnetic multilayers in current perpendicular to the plane of the 
layers geometry are defined by diffusive scattering in the bulk of the layers, and diffusive 
and ballistic scattering across the interfaces. Electrons scattered at the interfaces may 
be scattered again in the bulk of the layers, and, in order to better describe transport 
in the entire structure, the transport in the bulk of the layers and the transport through 
the interfaces should be treated self-consistently. Embedding the ballistic and diffusive 
interface scattering in the framework of the diffusive scattering in the bulk of the layers is 
the unifying idea behind the problems that I address in my work: the problem of finding 
the interface resistance and the problem of current-induced magnetization switching. 

1.4.1 Interface resistance 

Chapters |3] and m of this work are devoted to the problem of finding resistance due to the 
interfaces in metallic multilayered structures. The presence of the specular and diffuse 
scattering at interfaces leads to a finite voltage drop across them, ifTTl ^ [THl [121 and, 
hence, to a resistance. An origin of interface resistance, and the procedure of finding a 
resistance of an interface between two ballistic conductors is discussed in Chap.|2| In 
the realistic multilayered structure, where transport in the bulk of the layers is diffusive, 
resistance of the interfaces has to be incorporated into the resistance of the whole system. 
In order to do that, the resistors-in-series model was developed, and is currently used to 
analyze experimental data. lEUlETlE^E^ Within this model, resistance of an interface 
is treated independently of the resistance of the bulk of the layers, and total resistance 
of a multilayered system is a sum of the bulk and interface contributions .(Fig . 1 1 . 91) 

But the assumption that the resistance of the interface and resistance of the bulk are 
independent is not correct. An electron reflected from the interface can be scattered in 
the bulk and go back to the interface, not to the reservoir, hence changing the current 
through the interface. In terms of chemical potentials it means that in addition to an 
abrupt drop at an interface, there is a gradual change in the chemical potential within a 
mean-free-path from the interface. [24, 25, 26, 27, 28,29] Fig. ll.lOl shows schematically 
a chemical potential profile at an interface between two semi-infinite metallic layers. 
In Fig. II.IOL as well as in the all following pictures of chemical potentials, only the 
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Figure 1.9: Resistors-in-series model, pi and p2 are resistivities of the metallic layers 
Ml and M2, di and d2 are thicknesses of Mi and M2, Rint is a resistance of the interface 
between Mi and M2. 



potential drop due to the scattering from the interface is shown; the drop due to the 
resistivity of the adjacent layers is not shown. The resistance measured far from the 
interface, Ri, is different from the one that would be measured in its immediate vicinity, 
Rq. In the multilayers with layer thicknesses of the order of electron mean-free path, 
the chemical potential doesn't reach its asymptotic value within a layer. This leads to a 
breakdown of the resistors-in-series model. While the resistance of the whole system is 
still found by adding resistances due to the bulk of the layers and resistances due to the 
interfaces, interface resistances are not independent of the properties of the bulk of the 
layers, in particular, of the ratio of the layer thickness to the electron mean-free path in 
this layer. 

In order to determine the contribution of the interfaces to the resistance of the whole 
multilayered structure with layer thicknesses of the order or less the electron mean-free 
path, the form of chemical potential everywhere in the system has to be considered, 
and there was no such study for the systems with both specular and diffuse scattering 
at interfaces. In Ref. ll25ll the variation of the chemical potential was studied, however 
without examining how diffuse scattering at the interface altered its pattern, while in 
Ref. [|30| the interface resistance Ri (their Ra/b) coming from the specular reflection 
at the interface together with diffuse scattering in the bulk of the layers is calculated far 
from the interface. More recently the change of interface resistance Ra/b with diffuse 
scattering at the interface was calculated in Ref. ll3T1l . however only far from the inter- 
face. Contrary to the works ll25ll30ll3T1l . where ab-initio calculations of the interface 
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Figure 1.10: Chemical potential profile in two semi-infinite metallic layers. Only the 
potential drop due to the scattering from the interface is shown; the drop due to the 
resistivity of the adjacent layers is not shown. Rq is the resistance measured directly at 
the interface; Ri is the resistance measured far from the interface. 

resistance due to specular reflections at interfaces are performed, simple free-electron 
model for reflections at the interfaces is used in this work. In Chap.|3J a procedure allow- 
ing one to find a chemical potential profile everywhere in the multilayered structure in 
CPP geometry, including the close vicinity of the interfaces, and take into account both 
specular and diffuse scattering at the interfaces will be developed. (3T\ Results of the 
calculation of the resistances due to the interfaces in various systems will be presented 
in Chap.H 

1.4.2 Current-induced magnetization switching 

Chapters 13 and 1^1 of this work are devoted to the problem of the magnetization switching 
that is driven by spin-polarized current in noncoUinear magnetic multilayers. It is known 
that the transfer of spin angular momenta between current carriers and local moments 
occurs near the interface of magnetic layers when their moments are non-coUinear. The 
specular scattering of the current at interfaces between magnetic and nonmagnetic lay- 
ers that is attendant to ballistic transmission can create spin torque. iT, "Fl 191 [TDl [TTl 1^ 
154. .35J However, to determine the magnitude of the transfer, one should calculate the 
spin transport properties far beyond the interface regions. [Il36| In this work, I consider 
the effect that the diffuse scattering in the bulk of the magnetic layers and the diffuse 
scattering at interfaces have on the spin-torque; the spin transfer that occurs at interfaces 
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is self-consistently determined by embedding it in the globally diffusive transport cal- 
culations. ll37l The ballistic component of transport can be accommodated within the 
formalism described below, but this requires the knowledge of the band structure in the 
layers and is outside the scope of this study. ll38l 

A model system to calculate the spin torque is a magnetic multilayer whose essential 
elements consist of a thick magnetic layer, whose primary role is to polarize the current, 
a thin magnetic layer that is to be switched, a nonmagnetic spacer layer so that there is no 
interlay er exchange coupling between the thick and thin layers, and a nonmagnetic layer 
or lead on back of the thin magnetic layer; see Fig. ll.l l| The transport in the multilayer 
is considered as a diffusive process, and the interfaces are taken into account via the 
boundary conditions. In order to discuss the angular momentum transfer between the 
spin-polarized current and the magnetic background, one has to consider the exchange 
interaction between the accumulation and the background, or "sd" interaction, described 
by the Hamiltonian operator Hint = — Jm ■ M.^. The term {J /h)m x M^^ can be shown 
to exist in the equation of motion of the spin-accumulation by considering the quantum 
Boltzmann equation for the distribution function which takes the spin of an electron into 
account. [39 1 It describes a deterministic or ballistic precession of the accumulation due 
to the "sd" interaction when the magnetization directions of the spin-accumulation and 
the local moments are not parallel. The exchange coupling parameter J is the difference 
between the energies of the spin-up and spin-down electrons for each particular value of 
the electron momentum, J{k) = ei{k) — ei{k) (see Fig. ll.l2l ). averaged over the Fermi 
surface. The value of J can be obtained from the band structure calculations. For cobalt, 
I use the value of J =0.3 eV. I.40J For permalloy, J has been directly measured to be 
about 0.1 eV. MB 

The exchange interaction leads to the appearance of the component of the spin ac- 
cumulation transverse to the local magnetization direction. The transverse spin accu- 
mulation produces two effects simultaneously: one is to create a magnetic "effective 
field" acting on a local magnetization in the layer, and the other, the "spin torque", is to 
increase or decrease the angle between the magnetizations in the ferromagnetic layers. 
Contrary to the previous treatments, |I71[8l|9l[T0l|TT| both these effects enter the equation 
of motion for the local magnetization (the Landau-Lifshitz-Gilbert equation, Eq. (15. 81) ) 
on an equal footing. 

One of the consequences of the global approach to evaluate the effect the polarized 
current has on the background magnetization is that due to the presence of the long 
longitudinal spin diffusion length, longitudinal and transverse components (to the local 
magnetization direction) of spin accumulations are intertwined from one layer to the 
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Figure 1.11: Multilayered pillar-like structure used for current induced reversal of a 
magnetic layer. FM2 is a thick ferromagnetic layer with the thickness exceeding A^; 
and local magnetization M^^-* = cos ^e^ — sin 9ey, Sp is a thin nonmagnetic spacer, FMl 
is a thin ferromagnetic layer with the thickness tp and local magnetization M^^^ = e^, 
and NM is a nonmagnetic back layer. 
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Figure 1.12: Exchange splitting 
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next. As will be shown in Chap.|3 this leads to a large amplification of the spin torque 
acting on the free ferromagnetic layer. The angular momentum transferred to a thin 
layer /ar exceeds the transverse component (to the orientation of the magnetization of 
the free layer) of the bare portion of the incoming spin polarized current, i.e., the part 
proportional to the electric field. This, in turn, may lead to the reduction of the critical 
current necessary to switch spintronics devices (see Sec. II. 31) . 1 16| 

The transfer of angular momentum from the polarized current has an effect on the 
voltage drop across the multilayer being studied. |I3 The experimental data on several 
multilayered structures has confirmed that there are corrections to the simple cos^(6'/2) 
dependence of the CPP resistance, where 6 is the angle between the magnetization di- 
rections in the magnetic layers. L42>i43J The formalism described in this work allows to 
evaluate these corrections. 

In Chap. 121 1 will present a spin transfer model in which the equation of motion of the 
spin accumulation (essentially a spin-diffusion equation) is solved in order to describe 
the effect of the spin-polarized current on the background magnetization. Based on this 
equation, a simplified system of two thick magnetic layers separated by a non-magnetic 
spacer, and, in Chap. a realistic multilayered structure depicted in Fig. II. Ill will be 
studied. 
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Chapter 2 

Description of transport in multilayers 



In this chapter, the different formalisms for calculating electrical transport in solids will 
be reviewed. There are three most widely used ways to describe transport - semiclassical 
Boltzmann equation approach, Landauer formalism, and Kubo approach. 

The semi-classical Boltzmann equation approach combines Newton's law with a 
probabilistic description of random scattering forces. It can be used to describe transport 
in the systems where the quantum interference effects do not play a significant role, i.e., 
the electron mean-free path is large compared to the lattice constant, and small relative 
to the macroscopic dimensions of the system. The Boltzmann equation approach gives 
conductivity in terms of classical parameters of electrons - velocity, mean-free path, 
density. This formalism will be discussed in Sec. 12. II below in detail as it is the main 
tool used to obtain results in the present work. 

In the Landauer approach, ll44l l45ll current flow can be viewed as a consequence 
of the injection of carriers at contacts of a sample, and their probability of reaching 
the other end. In this approach, a current through a conductor, and conductivity, are 
expressed in terms of scattering properties of the system, or in terms of the probability 
that an electron can transmit through a conductor. P31 W7\ The Landauer formula for 
the conductance of a localized scatterer (an interface between two metallic layers, for 
example) will be derived and discussed below in Sec. 12.21 

The Kubo formalism is a method of calculating the response of a many-particle 
system to an external potential, for instance, the electrical current in response to an 
electric field. The external field is treated as a small perturbation on the equilibrium state 
of the system, eliciting a linear response, whose magnitude measures the corresponding 
transport coefficient, which is given in terms of the equilibrium properties of the system. 
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i.e., in zero field. The electrical conductivity tensor for example, may be expressed 
abstractly by the Kubo formula: [48 1 

1 

The conductivity depends on the time correlation between a component of the current 
operator ji/(0) at time zero and the component j)j(t) at some later time t, integrated 
over all time and evaluated as the average of the expectation value of the product over 
the equilibrium ensemble. The direct application of the Kubo formulas can be cum- 
bersome, but they provide an exact basis for various theorems involving the transport 
coefficients. Fortunately, all three approaches - Boltzmann equation, Kubo, and Lan- 
dauer (for electrical conductivity) - can be shown to be equivalent, and yield the same 
results for transport coefficients in some limits. I.47J 



2.1 Boltzmann equation 

The Boltzmann approach assumes the existence of a distribution function /(k, r), which 
is a local concentration of carriers in the state k in the neighborhood d^r of the point r 
in space. The total number of the carriers in the small volume d^rdk of the phase space 
is 

2 X {4^ of the allowed k states per unit volume in the momentum space) 
X dkx /(k, r)d^r, 

where the factor of 2 enters because of the spin degeneracy. The number of allowed 
states per unit volume dk is inverse of the volume occupied by one state in the momen- 
tum space. As follows from the Schrodinger equation for a free electron with periodic 
boundary conditions in the Lz box, the allowed k-states are spaced within 

2n / Li in each direction, so that the volume in the k-state per one state is {2txY /V . The 
total number of carriers in the phase volume (i'^rc/k is 

d T 

diV = 2/(k,r)^^ik. (2.1) 

The distribution function can change due to the diffusion of the carriers, under the 
influence of the external fields, and due to scattering. The Boltzmann equation states 
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that in a steady state at any point, and for any value of k, the net rate of change of 
/(k, r) is zero, i.e. 



The rate of the change of the distribution function due to diffusion (the motion of the 
carriers in and out of the region r) can be found as follows. Suppose that Vk is the 
velocity of a carrier in state k. Then, in an interval t, the carriers in this state move a 
distance v^t. Since the volume occupied by particles in phase space remains invariant 
(Liouville's theorem), the number of carriers in the neighborhood of r at time t is equal 
to the number of carriers in the neighborhood of r — Vk^ at time 0: 



This means that the rate of the change of the distribution function due to diffusion is 



External fields change the k-vector of each carrier, according to Newton's law, at the 
rate 




where F is the force acting on the carriers due to, for example, an electric field E, so 
that F = eE. According to Liouville's theorem in k-space, one can write 



so that the rate of the distribution function change due to the external fields is 




/(k,r,t) = /(k,r-vkt,0). 




(2.2) 



/(k,r,i)=/(k-kt,r,0), 




(2.3) 



The rate of change of f(k,r) due to scattering is 




[/(k',r)(l-/(k,r))Pk',k-/(k,r)(l-/(k',r))Pk,k']c?k'. (2.4) 
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The process of scattering from k to k' decreases /(k, r). The probability of this process 
is proportional to the number of carriers in the initial state k, /(k, r), and to the number 
of vacancies in the final state k', 1 — /(k', r). It is also proportional to the scattering 
probability Pk,k', which measures the rate of transition between the states k and k' if the 
state k is known to be occupied and the state k' is known to be empty. There is also the 
inverse process, the scattering from k' into k, which increases /(k, r). It's probability 
is proportional to /(k', r)(l — /(k, r)). The transition rate from k' into k is Pk',k- The 
summation over all possible k' states has to be performed. 

Combining the equations (12.21) . (12. 3I) . and (I2.4I) . one obtains the Boltzmann equation 
for the distribution function: 

Vk- Vr/(k,r) + |E- Vk/(k,r) = |[/(k', r)(l - /(k, r))Pk',k (2.5) 

- /(k,r)(l-/(k',r))Pk,k']rfk'. 

In this form, the Boltzmann equation is a nonlinear integrodifferential equation. In gen- 
eral, Pk,k' may depend on the distribution function /(k, r), and on the distribution of the 
scatterers. The nonlinearity may be removed provided that the principle of microscopic 
reversibility is valid, or the symmetry Pk,k' = Pk',k exist. This is usually the case if 
the crystal and scattering potentials are real and invariant under spatial inversion. [49| If 
the scatterers are sufficiently dilute and the potential describing the interaction between 
a carrier and a scatterer is sufficiently weak, Pk ^' is independent of the distribution 
function /(k, r). 

The Boltzmann equation may be simplified to a linear partial differential equation 
by the relaxation-time approximation. This approximation assumes that there exist a 
relaxation time r such that an electron experiences a collision in an infinitesimal time 
interval dt with probability dt/r. In general, the collision rate 1/r may depend on the 
position and the momentum of the electron: r = r(k, r). The additional assumptions 
are necessary to express the fact that collisions drive the electronic system toward local 
equilibrium. First, the distribution of electrons emerging from collisions at any time 
is assumed not to depend on the structure of the nonequilibrium distribution function 
/(r, k, t) just prior to collision. Second, if the electrons in a region about r have the 
equilibrium distribution function 

/(r, k, t) = r (r, 6) = , (2.6) 

g fcsT(r) _j_ I 
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where both the chemical potential fi and the temperature T may depend on the coordi- 
nate, the collisions will not alter the form of the distribution function. |50J Assuming 
a uniform temperature distribution in space, the equilibrium distribution function takes 
the following form: 

f{rXt) = f{e) = ^ , (2.7) 

e'^BT +1 

where ej? is an electron Fermi energy. 

It can be shown ll49l that in the relaxation-time approximation the collision term in 
the Boltzmann equation simplifies to 

/(k,r)-/°(6) 
r(k,r) • 

This form of the scattering term reflects the fact that the role of the collisions is to 
bring the system to an equilibrium. The Boltzmann equation in the relaxation-time 
approximation takes the form 

Vk ■ Vr/(k, r) + ■ Vk/(k, r) = .ll^lh^M, (2.8) 

The relaxation-time approximation provides the same description as the full Boltzmann 
equation if it is applied to spatially homogeneous disturbances in an isotropic metal with 
isotropic elastic scattering. |49| In this case, the relaxation time is defined as 




r(k) 

This relaxation time is called the transport relaxation time. 

Within the relaxation-time approximation, the Boltzmann equation (12.81) may be 
solved, and the electrical conductivity, which is the proportionality constant between 
the electrical current and the field, j = crE, may be found. In the homogeneous medium 
kept at constant temperature, the terms proportional to the spacial gradient V are equal 
to zero, and the distribution function /(r, k) turns out to be 

/(k) = r(e(k)) - -Lii^(k)vk ■ eE. (2.9) 
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Figure 2.1: Displacement of the Fermi surface in the momentum space under the influ- 
ence of the electric field 



This equation may be written as 

Assuming that the electron relaxation time does not depend on it's momentum, r(k) = 
r, it looks as if the whole Fermi surface had been shifted by the amount (er/^)E in k- 
space (Fig 12. II) . But, in fact, only the electrons close to the surface of the Fermi sphere 
have moved, the ones near the bottom of the conductivity band, deep within the Fermi 
sphere, are not really affected by the field. One can say that only the electrons with 
energies close to the Fermi energy, and, hence, the velocities close to the Fermi velocity 
participate in transport. The conductivity of a metal depends only on the properties of 
the electrons at the Fermi level, not on the properties of all electrons in the metal. 

The electrical current density j is defined as 



where /(v, r) is the density of the carriers in the (r, v) space, so that the number of 
electrons in the volume (Prdv is 



Comparing this expression with the equation (12.11) . the following relationship between 




(2.10) 
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/(v, r) and /(k, r) may be obtained: 



where m is the electron mass, and the relation mv = hk is used. Considering the case 
where both j and E are in the a;-direction, the relaxation time r is independent of the 
momentum, and substituting /(k, r) from the equation ( 12.91) . one obtains the following 
expression for the conductivity: 

1 £!lf!n-)7f_M£)Vw,, (2.12) 



where the fact that / Vk/°(e(k))(iv = is used, since /°(e(k)) is isotropic in k. The 
factor 1/3 appears in the expression since has been averaged over the Fermi surface, 
so that t>^ = |t>^. In a metal for the temperature close to 0°K, the function (— (9/°(e) / de) 
behaves like a delta-function at the Fermi level, 5{e — ep), and the expression for the 
conductivity (12.121) reduces to 



ki e T ne T 



m 



(2.13) 



where n = kp/Sn'^ is the electron density, and m is the electron mass. Strictly speak- 
ing, the above expression for the electrical conductivity is only valid for free electrons. 
The free-electron description of a metal will be used throughout this work, so the equa- 
tion (12.131 ) will be used to calculate the conductivity. 

An alternative form of the scattering term in the Boltzmann equation may be ob- 
tained as follows. Assuming the rate of transition between the states k and k' to be 
independent of k and k', Pk,k' = -Po, the scattering term Eq. (12.41) may be written as 

fdf\ /(r,k)-/(r) 



dt I ., T 

/ scatt 



where 



1 

r 



Podk, 



m = ^f f{r,k)dk. (2.14) 

4:71 Jfs 
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The linearized Boltzmann equation then takes the form 

,„,.Ef!..,.V/(r.k)^-/(^J^l^. (2.15) 

oe T 

Introducing g{r, k) such that /(r, k) = /°(e(k)) — ^g(r, k), the function / (r) may be 
written as 

de 4:11 Jfs 

and the equation for the function g{r, k) takes the form 



/(r) = r(e(k)) - ^- / ^(r,k)dk, (2.16) 

oe Air Jfs 



■ Vg{r, k) + = ev^-E+^, (2.17) 

r r 



where /i(r) is defined as 



/i(r) = — / r/(r,k)rfk, (2.18) 

Air Jfs 



so that g(r, k) may be written as 

^(r,k) = /i(r) +^(r,k), 

where 

/ ^(r,k)rfk = 0. (2.19) 

Jfs 

The function fi{r) defined in Eq. (12.181) is the non-equilibrium part of the chemical 
potential, since the distribution function /(r, k) may be written as 

/(r, k) = r (6(k)) - ^/.(r, k) - -I-~g{r, k) ^ -j-~g{r, k) 

(2.20) 

Equations (12.171) and (12.181) are the starting point of Chap.|3l where the equations for the 
chemical potential profile in the multilayered metallic structures are derived. 

To conclude the discussion of the Boltzmann approach to the electrical transport in 
solids, the diffusion equation on the concentration of carriers will be derived below. The 
non-stationary linearized Boltzmann equation takes the following form: 

df{x,k,t) , df{x,k,t) dfM /(x,k,t)-/» 

^ ^ e E^v^ = , (2.21) 

ot ox oe T 
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where, for simplicity, electrical field is assumed to have only one component, E^, and 
spacial gradient in s-direction only is considered. Integrating Eq. (I2.21I) over all v- 
space, and using the definitions of the concentration and current density 

n{r,t)= J /(r,v,t)rfv, (2.22) 

and 

i = el v/(r,v,t)rfv, (2.23) 
where m is the electron mass, the continuity equation 

dn(x,t) IdiJx.t) ^ 

— ^ \ = (2-24) 

dt e ox 

may be obtained. The term proportional to / Vx{df /de)dv can be shown to be zero, and 
the term proportional to /(/ — /)(iv is zero due to the condition (12.191) . Multiplying the 
Boltzmann equation ( 12.211 ) by and integrating over all v-space again, the following 
equation is obtained: 

e dt J ^ dx 4n^\hJ J ^ de 

= (2.25) 

e r 

The term proportional to / Vxf^dv can be shown to be zero. Assuming the electrical 
current to be independent of time, and writing the second term, / vl(df/dx)dv, as 
< vl > J fdv = Ivp-^ J fdw, the following expression for the current is obtained: 

Ux)=aE,~eD^^, (2.26) 

where D = \vpT is the diffusion constant, a = is the electrical conductivity. 

Equation (12.261) states that within a linear response model, the current is proportional 
to the electrical field and the gradient of the carriers concentration with the negative 
sign. Substituting the equation (12.261) into the continuity equation (12.241) . the diffusion 
equation is obtained: 

dn{x,t) ^ ^ d'^n{x,t) 
dt dx"^ 
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Equations (12.241 ) and (12.261 ). generalized to the spinor form, are the starting point of 
Chap.lH where the equations for the spin-accumulation are obtained. 



2.2 Landauer approach - interface resistance. 

The conductance G due to elastic scattering of an obstacle, characterized by transmis- 
sion and reflection coefficients T and R, is given by 

T 

G = — -, (2.28) 

vrn. R 

where e is electron charge, and h is Plank's constant. At zero temperature T and R are 
evaluated at the Fermi energy. Equation (12.281) applies to samples of arbitrary shape and 
structural complexity, and may be used, for example, to calculate the conductance of a 
planar barrier, such as an interface between two metals. Below, equation (12.281) will be 
derived within an approximation of a single conductivity channel, PtI 151114611521 which 
makes a system effectively one-dimensional. Expression (12.281) can be generalized for 
the case of many independent conducting channels, and for nonzero temperature, iBTl 
Wl\ but these generalizations won't be covered here. 

Consider a barrier connected through ideal ID wires to some external source (a 
pair of reservoirs with different chemical potentials jii and 112) which drives a current 
/ through the system. The barrier is characterized by its transition coefficient T and 
reflection coefficient R = 1 — T . (Fig. 12.21 ) 

Consider first two reservoirs without a barrier. Left reservoir emits electrons with 
energies up to a quasi-Fermi energy (chemical potential) fii. Right reservoir emits 
electrons with energies E up to the chemical potential /i2. Energies are measured relative 
to equilibrium chemical potential (common for both reservoirs). Current emitted by the 
left reservoir and going to the right is 

= —evniE — /ii) —evno(E) + ev— — fii, 

oE 

where n{E) is electron density per unit length, v is the velocity component along the 
wire at the Fermi surface. Current emitted by the right reservoir to the left is 

J- = evn{E - ^2) ~ evnQ{E) - 

oE 
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Figure 2.2: One-dimensional barrier with the transition coefficient T and reflection co- 
efficient i? = 1 — T, connected to an external source (a pair of reservoirs with different 
chemical potentials /ii and /i2.) 



Then, total current is given by the sum of and / and is equal to 

I = ev—{fii-fi2)- 

dn/dE is the density of states for two spin directions and for carriers with positive 
velocity. In one-dimensional case, dn/dE = l/nhv ll53l and 

I = -^(/ii -/^2). 

This is the total current emitted by the left reservoir due to the difference in the quasi- 
Fermi levels. Carriers have a probability T for traversal of the sample and a probability 
R of being reflected back. Assuming that T and R don't depend on the electron energy 
in the energy range fxi > E > ji2, the net current flow may be written as 

/ = -^T(/ii -/ia). 
irn 

The difference in chemical potentials of the two reservoirs yUi — /i2 can be identified with 
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the drop of electron potential eV across the source. Conductance, given by G = I /V is 



G ^ T 

(/ii-/i2)/e nh 

and resistance of the system is 

= (2.29) 

Expression ( I2.29I ) gives a non-zero resistance in the case of completely transparent bar- 
rier, T = 1. The resistance in a system consisting only of the ideal wires and reservoirs 
is a finite quantity 

■jrfy 

Gc^ = ^ = 12.9 kn. (2.30) 

This quantity is called the contact resistance. It occurs between a large reservoir with 
the infinite number of states that an electron can have, and a narrow wire with only 
one conductivity channel in the same way as traffic jam occurs when 4-lane highway 
transforms into 2-lane street. If a larger number M of the conductivity channels is 
allowed to exist in a wire (when the wire is three-dimensional but still very thin, so 
that the perpendicular motion of an electron is quantized), the contact resistance can be 
shown to have the form [47 1 

e2 M' 

The wider the wire, the larger the number of transverse modes M can propagate in it 
and the smaller the contact resistance is, and it is negligible for thick contacts. 

The resistance of the barrier alone can be calculated by subtracting the contact re- 
sistance (equation (12.301) ) from the total resistance of the "barrier-contacts-reservoirs" 
system (equation (12.291) ). and is equal to 

G" = (2.31) 



this is equivalent to the expression (12.21) for the conductivity when one uses 1 — T = R. 

Resistance is usually associated with energy loss. In the considered model, scattering 
both at the barrier and in the wires is perfectly elastic {R + T = 1), so that there seems 
to be no energy loss in this system. Nevertheless, resistance exists. The resistance of the 
barrier can't be considered without connecting the barrier to the reservoirs. Electrons 
reflected from the barrier go back to the reservoir. Instead of reissuing the electron that 
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entered the reservoir, a new one comes out, so that information about the momentum of 
the electron is lost in a large reservoir, the outcoming electron does not have the same 
direction of the momentum as the incoming one had. This loss of information is the 
origin of the resistance of the barrier. ll54l The existence of resistance in the Landauer 
approach requires the presence of the reservoirs, but its magnitude depends only on the 
elastic events at the barrier. 
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Chapter 3 

Interface resistance in multilayers - 
theory 

The Landauer formula ( 12 .291) is derived for an obstacle connected to the source via ideal 
scattering-free wires. The purpose of this work, though, is to find a resistance due to 
interfaces in a realistic multilayered structure consisting of several diffusive metallic 
layers and interfaces between them. As discussed in Chap. [H scattering at the inter- 
faces is affected by the presence of the diffuse scattering in the bulk of the layers, and 
contributions from the bulk of the layers and the interfaces to the total resistance of the 
structure can't simply be added. In this section, the procedure allowing one to find a 
chemical potential profile and resistance due to the interfaces in a multilayered metal- 
lic structure in CPP geometry that takes into account an interaction of the scattering at 
the interfaces and scattering in the bulk of that layers will be developed. Transport in 
the bulk of the layers will be described by semiclassical Boltzmann equation, and in- 
terfaces will be taken into account via boundary conditions. Both specular and diffuse 
scattering at the interfaces will be considered. While Landauer's result for the interface 
resistance won't be explicitly used in this work, only scattering properties of the inter- 
faces, their reflection and transmission coefficients will be used, which is a signature of 
the Landauer approach. 

First, a general solution of the semiclassical linearized Boltzmann equation in CPP 
geometry is presented. Next, boundary conditions at an interface between two layers in 
the presence of both specular and diffuse scattering at the interface will be discussed. 
Then, a system of integral equations describing a chemical potential profile in two semi- 
infinite metallic layers will be derived. 0211 Integral equations on the chemical poten- 
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tials in the systems consisting of three and five layers are presented in Appendices IA.4I 
and IA.5I This chapter is completed with determining the special form of diffuse scat- 
tering such that resistances measured at the interface and far from it are the same for 
a system consisting of two semi-infinite layers of identical metals and an interface be- 
tween them. The forms of diffuse scattering for which a prediction can be made as to 
whether the potential drop at the interface is bigger or smaller than that measured far 
from it will be indicated. 



3.1 General solution of the Boltzmann equation 

In the relaxation-time approximation the linearized Boltzmann equation in CPP geom- 
etry takes the form (see Eq. (12.171) and ji(z) denoting the non-equilibrium part of the 
chemical potential): 

dg{e,z) g{e,z) , K^) 

— 7; \ = e-&2 H , (i.l) 

OZ VzT VzT 

where e is an electron charge, r is an electron relaxation rate, and are an electron 
velocity and an electric field in the direction z perpendicular to the plane of the layers, 
and 9 is the angle that the electron momentum k makes with z axis which takes values 
from to TT. Chemical potential ii{z) is the average of g{9, z) in momentum space (see 
Eq. (ITTBb V 

^(z) = - r g{9, z) sin 9d9. (3.2) 
2 Jo 

Equation (13.11 ) is an ordinary first-order differential equation. Since the chemical 
potential jj,{z) enters the right hand side of this equation, its solution (7(6*, z) in the region 
between zi and Z2 of z axis will depend on via an integral: 

g{9, z) = r dz'k{9, z, z')fi{z') + y{9, z), (3.3) 

where k{9, z, z') and y{9, z) are functions defined by the geometry of the problem and 
the boundary conditions. By using the definition of the chemical potential (13. 21) . an 
equation for can be obtained by integrating equation (13.31) over 9: 

fi{z)= dz'K{z,z')fi{z') +Y{z), (3.4) 

J Z\ 

where K{z, z') = \ k{9, z, z') sin 9d9, Y{z) = \ £ y{9, z) sin 9d9.^ 
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Figure 3.1: Scattering at an interface: a) schematic representation of electron scattering, 
b) direction of electron fluxes at an interface. 

In the structure consisting of several layers (N), boundary conditions connect distri- 
bution function in a layer (Lj) with the distribution functions in the adjacent layers (Lj). 
Equation (13.41) can then be generalized to describe the chemical potential everywhere in 
the system: 



The integral equation (13.51) is a system of Fredholm equations of the second kind, 
which can be numerically solved to obtain the chemical potential profile everywhere 
in the layered structure. [.32.. .561 The functions Kij{z, z') and Yi{z) are defined by the 
boundary conditions at the interfaces. Numerical procedure of solving the Fredholm 
equation of the second kind is presented in the Appendix lA. 61 

3.2 Boundary conditions at the interface between two 
layers 

At the interface between two layers, an incoming electron is scattered via reflection and 
transmission. There is also diffuse scattering in all directions. In the present work, a 
simple model for the diffuse scattering will be used, where a single parameter S rep- 
resents an amount of electrons which are not scattered diffusely, so that the specularly 
reflected and transmitted electron fluxes are scaled by S*. 1 — S* is the amount of electrons 
scattered diffusely. Figure l3TTk ) represents schematically how an electron is scattered at 
an interface. 




(3.5) 
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The boundary conditions for the distribution function / at the interface 2 = be- 
tween two semi-infinite layers in the presence of diffuse and specular scattering take the 
following form: ||24ll57t 



S{vi, V2)R2l{Vi, V2)f{v^, 0) + S{Vi,V2)Ti2{Vi, V2) f {v^ , 0) 

+ [l-Siv,,V2)]F{f), 



< 



(3.6) 



S{vi, V2)Ri2{vi, V2)f{v^, 0) + S{vi,V2)T2i{vi, V2) f {vf , 0) 
+ [l-S{v„V2)]FU). 



where 



F{f) = -( c/k[K|/(t;<,0) + \v>^f{v>M. 

\L JfS 



(3.7) 



n = [ dk6{e-eF)\v, 



JFS 



ep is the electron Fermi energy, vi {V2) is the electron velocity in the layer 1(2), and 
{v"^) refer to f ^ > {vz < 0). Rij{vi,V2) (Tij{vi,V2)) is the reflection (transmission) 
coefficient at an interface for an electron travelling from the layer i to the layer j,and 
1 — S{vi,V2) describes the diffuse scattering. As shown at the Fig. 13.1b ). an outgo- 
ing particle flux /(u^) (/(wf )) comes from the specularly reflected particle flux (first 
term), the specularly transmitted particle flux (second term), and the particle flux from 
all directions diffusely scattered in the direction of V2 (f f ) (last term). 

Similar boundary conditions, but without the last term, were investigated by Hood 
and Falicov for currents in plane of the layers geometry. iBTI In CPP geometry, cur- 
rent conservation requires the presence of the term responsible for the diffuse scattering 
at an interface in the boundary conditions. The proof of the fact that the boundary 
conditions (I3.6I) - (I3.7I) indeed conserve the current across an interface is presented in 
Appendix lA. 21 

At the interface between two metals, reflection and transmission coefficients Rij 
and Tij have to be found from quantum mechanical considerations, as reflection and 
transmission coefficients for the quantum mechanical flux. In general, Rij, Tij and S 
depend on both variables vi and t>2. In the quasi-one-dimensional layered structures, 
where the component of electron velocity parallel to an interface in conserved, only z 
component of velocity enters reflection and transmission coefficients. Furthermore, in 
the case of specular scattering, angle of reflection and angle of transmission are uniquely 
defined by the angle of incidence, so that Rij and Tij are the functions of the cosine of an 



35 



angle that an electron velocity makes with an interface only, Rij {cos Oi) and Tj^ (cos Oi), 
where 6i is the angle of incidence of an electron at the interface in i-th layer. Specular 
scattering at an interface between two metals is considered in detail in Appendix IA.3I 
In order to simplify the discussion of the effects of the diffuse scattering at an interface 
on the electron transport in multilayers, diffuse scattering parameter S is assumed to be 
constant, independent of the electron direction, with one exception, in Sec. 13.41 where 
the special forms of the diffuse scattering coefficient for which a shape of the chemical 
potential profile can be predicted analytically are considered. With these simplifications, 
the boundary conditions on the function g take the following form: 

r g{v>, 0) = 5i?2i(^^2)^?«, 0) + ST^2{eMv>, 0) + (1 - S)F{g), 

(3.8) 

[ g{v<, 0) = SRu{ei)g{v>, 0) + ST2i{e2)g{v<, 0) + (1 - S)F{g), 
where F{g) can be shown ( Appendix lA. 1.3 1) to have the form: 

2ii2 /•7r/2 

F{g) = — / ^K>,0)cos^sin^d^ (3.9) 

Vp2 + Vpi JO 

+ — / g(v^,0)cosesmede, 

where vpi and vf2 are the electron Fermi velocities in the first and the second layer 
correspondingly. 



3.3 Equations for chemical potential profile in a multi- 
layer 

In this section, solution of the Boltzmann equation in the system consisting of two semi- 
infinite metallic layers and interface between them (Fig. 13. 21) is presented, and equations 
on the chemical potential profile are derived. For simplicity, the same electron relax- 
ation time T is assumed in both layers. Solution of the Boltzmann equation (13.11) for 
corrections to the distribution function of the electrons moving to the right g{vi, z) and 
to the left f , z) in the first layer, and electrons moving to the right (?(f^, z) and to 
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Figure 3.2: System consisting of two layers. 



the left z) in the second layer take the following form: 



9{v2,z) 



^ ( K^) + /o dz' exp {eE., - ^ 

-/rrfz'exp(^) (eE 



\v,2\rl' 



(3.10) 



where yUi(2;) (/i2(-2)) is the chemical potential in the left (right) layer, Ezi (-E22) is an 
electric field in the left (right) layer, and Vzi {Vz2) is a component of an electron velocity 
in the direction z in the left (right) layer. Since only the z component of an electric field 
is present in CPP geometry, the subscript z will be omitted from E^i and Ez2 everywhere 
below. A, B are constants of integration, independent of z. They are defined by the 
boundary conditions (13. St . and equal to 



A = SRi2{ei)eEih\cosei\- ST2i{e2)eE2l2\cose2\ 
+ SRuiei) exp — —d^ 

J -00 \COS 61 J COS 61 



(3.11) 
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+ ST,,{e,) r exp ^di' + (1 - S)F{g) 

Jo \COS Ui J COS &2 

B = -SR2i{92)eE2l2\cos92\ + STi2{9i)eEih\cos9i\ 

+ SR2i{92) rexp(^]^de (3.12) 

JO \ cos 92 J COS 92 

+ ST^2{9i)f ex^(^\^di' + {l-S)F{g), 
j-oo ycos6'2 / costal 

where /j = vpiT is the electron mean free path in i-th layer, 9i is an angle that an 
electron velocity makes with z axis in i-ih. layer, and a dimensionless variable ^ = z/li, 
where i = 1 or 2 depending on in which layer the chemical potential is considered, 
is introduced. Since the boundary conditions connect the distribution functions at both 
sides of the interface, both angles 9i and 6*2 enter the expressions for A and B. From the 
definition of the chemical potential ( 13.21 ) it follows (see Appendix I A. 1.21 ) that 

2^i{z) = jQ^'^{g{vf ,z) + g{v<,z)) sin 9id9i 

(3.13) 

2/i2(2;) = Jo ^"^ {g{v^, z) + g{v^,z)) sin 92d92, 

where the averaging is performed over the electron momentum direction 9i and 6*2 in the 
first or the second layer correspondingly. Hence, the constants of integration A has to 
be expressed in terms of the angle 9i only, A(9i), and B in terms of the angle 92 only, 
B{92)- It can be shown (see Appendix lA. 31) that angles 9i and 92 are related; 9i can be 
expressed as a function of 92, 6*1(^^2), and vice versa, ^^2(^1)- It follows from this relation 
that Ti2{9i) = T2i{92) and Ri2{9i) = R2i{92). Note that when calculating the function 
F{g), change of angles is not required, which reflect the fact that in the process of diffuse 
scattering angles of incidence, reflection and transmission are not related, so a variable 9 
will enter the expression for F{g) without the subscript. The electrical fields Ei and E2 
in the layers are related due to the current conservation, so that aiEi = const throughout 
the system, where (jj is the conductivity of the layer i. Since electron relaxation time r 
is assumed to be the same in all layers, and the conductivity is inversely proportional to 
Vp, ll5^ electrical fields Ei and E2 are related via Vp^EJi = Vp2^2l2- 

The following equations on the chemical potentials in the left and right layers are 
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obtained from equations (I3.10I )- (I3.13I ): 



2/ii(0 = eE,h r 
Jo 



■/2 



X 



(1 + SRuiOi)) COS^i - 4^5Ti2(^l) COS02(^l) 

^F2 



X exp f — \ sinOidOi 
\ costal / 



J — oo 
r7T/2 

Jo 



^^P a — + S^mOi) exp — 

y cos Oi J \ cos Oi J 

■7r/2 



sm / exp 

\cosU J Jo \cosu 



X / exp 
Jo 



2(1 -g) r r^^'U^' 



|-7r/2 

X / exp 

JO 



(—^] sinOde f ' exp ( ^— 
\cosu J Jo \COSUj 



sin ^c^0 



(3.14) 



2/^2(0 



X 



-ei?2^2 



7r/2 



(1 + ^i?2l(^2)) COS ^2 - ^^^21 (^2) COS ^1(^2) 



X exp 
+ 







\ cos 612 / 



tt/2 







exp 



sin ^2C^^2 



cos 6*2 

7r/2 



+ 5i?2i(^2)exp 



COS ^2 



exp 



tan ^2 (^^2 

\ sin02 



cos 6*1 (6'2) / \ cos 6*2 / cos 6'i (6'2) 
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+ 



X 



X 



2(1 -S) 



1 + v'p2/'Vpi J -co 
7r/2 / ^' 



exp , 

\ COS 6 



2(1-5) 



r7r/2 / ^ \ 

sin 6d6 / exp sin 6d6 

Jo \ cos 9 I 



1 + ^Fl/^F2 -^0 

exp sin 6d6 / exp | | sin 6'ci6' 







cos 9 







cos^^ 



(3.15) 



These equations describe the chemical potential profile everywhere in the system 
consisting of two semi-infinite metallic layers in the presence of both specular and dif- 
fuse scattering at the interface between them. In a similar fashion, equations describing 
the chemical potential profiles in three- and five-layered systems can be obtained. Be- 
cause of the tedious calculations required for the derivation, only the final equations are 
presented in Appendices IA.4I and IA.5I Even in the simplest structure consisting of two 
layers of identical metals divided by an interface, the form of the chemical potential 
can't be found analytically, though is this case a certain prediction can be made. Results 
of the numerical solution of the equations on the chemical potential in various systems 
are presented in Chap. IH 



3.4 Special forms of diffuse scattering 

In the simple case of identical metals at both sides of an interface, special forms of dif- 
fuse scattering ^(cos 9) can be found, for which the chemical potential remains constant 
outside the interfacial region, so that the resistance measured at the interface is the same 
as the resistance measured far from it. In the identical layers, vpi = vf2, di = 02, and 
chemical potential is antisymmetrical around z = 0: 

/i2(0 = (3.16) 

Then, equation [3.14l for the chemical potential in the left layer, /Ui(^) for example, takes 
the form 

2/ii(0 = eE,k r^\l-S{9i) + 2S{9i)Ri2{9i))exp(-^]sin9icos9id9 
Jo \cos9i J 
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+ r fiiiOdC (3.17) 



oo 



X [exp ( I + 8(9,) (2i?i2(^i) - 1) exp ( ) ) tan^it/^i. 

JO \ \ cos oi / \ COS VI / / 

The derivative of the chemical potential takes the form 

2/i;(0 = eE^h r^\i-s{e,) + 2S{ei)Ru{ei))exp(^-]sme,de, 

Jo N^coso'iy 
° /ii(r)c?r (3.18) 

"OO 

X (2i?i2(^i) - 1) exp (^-±pi '-^de, 

Jo y costal J costal 

- / f^iiUd^ exp — — d^i 

J~oo JO \ costal / costal 

+ / ^^l[i)di / exp — —dOi. 

J^ Jo Y costal J costal 

The condition fi\ (^) = 0, so that the chemical potential in the left layer /ii ) is constant, 
can be satisfied when the reflection coefficient RuiO) and diffuse scattering parameter 
S{6) are related via 

^'^^^ = {cose + m-2R{e)y ^^-^^^ 

where /3 is a constant. Then the chemical potential in the left layer is equal to Hi{z) = 
(5eEili = const, and in the right layer iJ.2iz) = —PeE2l2 = const, where electrical 
fields El and E2, and electron mean-free paths li and I2 are the same in both layers. The 
sheet resistance across the interface ARq and measured far from it ARi (Fig. ll.lOt are 

where m is an electron mass, is an electron Fermi velocity, same in both layers. 

The special form of diffuse scattering (Eq. (13.191 )) which makes fx{z) constant for 
z ^ can be found for any reflection coefficient. In particular, if an interface between 
metals is modeled as a sheet delta function potential, so that the reflection coefficient 
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Figure 3.3: Diffuse scattering S{9). Region A: potential drop at the interface is bigger 
then potential drop far from the interface. Region B: potential drop at the interface is 
smaller then potential drop far from the interface. Region C: region of special forms of 
diffuse scattering. 

takes the form R{9) = a/{a + cos^O), ||58l it exists for any /3 from 1 to a. These special 
forms of S occupy region C in Fig. 13.31 If the function S{9) stays above this region 
for all 6 (region A in Fig. 13.31) then the resistance across the interface Rq is larger than 
resistance Rs measured far from it. When S{9) stays below this region for all 9 (region 
B in Fig. 13.31 ) then Rq is smaller then Ri. For S{9) occupying both regions C and A or 
(and) B, no prediction for the chemical potential shape can be made. 
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Chapter 4 

Interface resistance in multilayers - 
results 



In this chapter, the numerical results for the chemical potential profile and resistance 
due to the interfaces in the multilayered systems are presented. 

The systems consisting of two, three, or five metallic layers, i.e., with 1, 2, or 4 
interfaces, in which I consider electrons move across a step-like potential. Energies are 
referenced to the Fermi level so that ep = 0, and the effective mass of the electron 
is assumed to be independent of the material and spin-orientation. In the i-th layer, 
electrons move in a constant potential Vi = —mvfp/2, where Vip is the electron Fermi 
velocity in the i-th layer. The same electron relaxation time r is assumed in all layers, 
and an electron mean-free path k = vpiT is introduced. Current conservation requires 
electric fields Ei in the layers to be related by Vpj^Eik = const (see Sec. 13.31) . The 
reflection coefficient at the interface between two metals takes the following form (see 
Appendix lA. 3 1 ): 



Besides the specular scattering at the interfaces, diffuse scattering is present, such that 
the diffuse scattering parameter S (see Chap. 13) is assumed to be independent of the 
angle of incidence, and the same for all interfaces. 

Equations (l3TT4l) - (l3TT5l) . (lAl5l) - (lAl9t and (IX20h - (IX28l) yield the chemical poten- 
tials in each layer of the multilayered system, normalized to the voltage drop within a 




(4.1) 
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mean-free path in this layer ^i{z)/eEili. Resistance due to the interfaces is proportional 
to the drop of the chemical potential: 




(4.2) 



where AR is the sheet resistance due to an interface, e is the electron charge and j is the 
electric current. In terms of the chemical potential drop normalized to the voltage drop 
within a mean-free path in the layer An/eEl, the interface resistance may be written as 



The factor Gn'^h/e^kp, different for each material and each spin-orientation, has to be 
taken into account when deducing the resistance due to the interfaces from the chemical 
potential drops obtained by solving the equations for the chemical potential profiles. 

i?Af, where N is the number of interfaces in the system, = 1, 2, 4, denotes the 
resistance due to the interfaces, measured far from the interfaces: 



For the two-layered systems, resistance Rq measured directly at the interface at 2; = 
(Fig. ll.lOb is introduced: 



In this chapter, the systems consisting of two semi-infinite metallic layers are con- 
sidered. The influence of the diffuse scattering at the interface on the chemical potential 
profile, and the influence of the diffuse scattering in the bulk of the layers on the in- 
terface resistance are investigated. The breakdown of the resistors-in-series model is 
anticipated. Next, the prediction about the resistors-in-series model breakdown is tested 
in the systems consisting of three and five layers. The interface resistance dependence 
on the layers thicknesses relative to the electron mean-free path is obtained for different 
parameters S of diffuse scattering at the interfaces. Finally, realistic multilayers con- 
sisting of two, three and five layers of ferromagnetic and non-magnetic metals (Co-Cu 
and Fe-Cr systems) are considered. Resistances due to the interfaces Rn are obtained 
as a function of the layer thicknesses and the amount of the diffuse scattering at the 
interfaces. I investigate the deviation of the resistance -Rat in the three- and five-layered 
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Afi A/i / Afi Q-n'^h 
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systems with the layer thicknesses less than the electron mean-free path in the layers 
from the resistance NRi that such systems would have if the interfaces were far enough 
from each other so that the scattering from the neighboring interfaces would not inter- 
fere. 

The numerical procedure of solving the integral equations for the chemical potential 
profile dTUll - dTTSll . (lAlSh - dAlOh and (IX20b - (IX28b is discussed in the appendixlA!6l 

4.1 Two layers 

4.1.1 Same metals 

Equations (13.161) and (13.171) are solved numerically ll32l to obtain the chemical po- 
tential profile in the system consisting of two identical metallic layers separated by 
a delta function potential, so that the reflection coefficient takes the form R{9) = 
a/{a + cos^ 9). |58| Fig. 14. lb ) shows how the chemical potential profile changes with 
the amount of diffuse scattering when S* is a constant between (only diffuse scat- 
tering is present at an interface) and 1 (only specular scattering is present). Interface 
resistance grows as the amount of diffuse scattering decreases, [24J and the chemical 
potential varies within a mean-free path of the boundary. The resistance Ri measured 
far from the interface can be larger or smaller then that measured at the interface, Rq, 
(see Fig. 11.101) . depending on the value of S. Fig. 14.1b ) shows the difference between 
Rq and Ri relative to Ri. This difference can be as large as 18-19%. When there is no 
diffuse scattering at the interface, S = 1, the difference between the resistance Ri and 
Rq, agrees with the results obtained by Penn and Stiles \.26l . 

4.1.2 Different metals 

Next, a system consisting of two different metallic layers and an interface between them 
is considered. [ |32| Figure l4!2l shows the step-like potential experienced by electrons in 
such a system. For z < 0, electrons move in a constant potential Vi = —mvpi/2, and 
V2 = -mvp2/2) for 2; > 0. 

Equations (13.141 ) and (13.151 ) with the reflection coefficients (14.11 ) are solved in order 
to investigate how the chemical potential profile and interface resistance change as the 
functions of the potential barrier height V2/V1, and the amount of diffuse scattering at the 
interface. In Fig. 14.31 the forms of the chemical potential for V2/V1 = 2, and different 
values of the diffuse scattering parameter S are shown. As in the case of the same metal 
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Figure 4. 1 : Chemical potential and interface resistance in two layers of identical metals: 
a) chemical potential profile for 5" = 0, 0.1, . . ., 1, normalized with respect to potential 
drop eEl within a mean-free path I of the metal, b) the difference between Rq and Ri 
relative to Ri. 
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Figure 4.2: Step-like potential experienced by electrons at an interface between two 
metals. 
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Figure 4.3: Chemical potential profile in two-layered system for different values of 
diffuse scattering S* = 0, 0.1, . . ., 1 at the interface, normalized with respect to potential 
drop eEili within a mean-free path k of each metal, so that the plots are anti-symmetrical 
around z = 0. 



on both sides of the boundary, the chemical potential varies within several mean-free 
paths of the interface. Resistance measured far from the interface, Ri, is bigger then that 
measured right at the interface, Rq, for interfaces where diffuse scattering dominates 
(S is close to 0). Ri is smaller then _Ro for the interfaces where specular scattering 
dominates (S is close to 1). Unlike the case of two identical layers, the resistance of the 
step-like barrier decreases as the amount of diffuse scattering decreases. [24| 

Figure B^ demonstrates the effect of diffuse scattering in the bulk of the layers on the 
interface resistance. In the case of a completely specular interface, 5=1, the variation 
of interface resistance as the function of the barrier height V2/V1 is shown in comparison 
with that obtained by Barnas and Fert in the system of ballistic layers. lITSl [T9| In the 
presence of diffuse scattering in the bulk of the layers, both Rq and -Ri are less than 
the interface resistance obtained for purely ballistic transport. In the diffusive layer, 
electrons reflected off the barrier can be scattered back on the barrier, and again have a 
chance to be transmitted through, making the interface resistance smaller. 

Figure 14.51 summarizes the behavior of the interface resistance as the amount of 
diffuse scattering at the interface, and the height of the potential barrier experienced by 
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Figure 4.4: Effect of the diffuse scattering in the bulk of the layers on the interface 
resistance in two-layered system. Red line is the Bamas and Fert result for the sheet 
interface resistance, ARbf, adapted from Refs. [il8l fT9^. for the ballistic transport in 
the bulk of the layers and specular scattering at the interface; blue and green lines are 
resistances ARq and ARi for diffuse transport in the bulk of the layers and specular 
scattering at the interface, S = 1, as the functions of the electron potentials ratio a) 
V2/Vi < l,h)V2/Vi > 1. 
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the electrons at the interface change. Arbitrary units for the resistance are used, since 
the calculations are performed irrespective of the particular metals. Note that the ratio 
of the electron Fermi velocities in the layers vf2/vfi is used as a variable, not the ratio 
of the electron potentials V2/V1, as in the Fig. 14.41 In Figs.|43^)-b), the variation of the 
interface resistance as the function of the amount of diffuse scattering at the interface 
for the system with vf2/vfi = 2 is shown. The resistance Rq measured directly at 
the interface is larger than the resistance Ri measured far from it, and Rq grows as the 
interface becomes less diffusive {S approaches 1), while Ri decreases slightly, so that 
the absolute value of the difference _Ri — _Ro grows, reaching about 70% if calculated 
relative to i?o, and 120% if calculated relative to Ri for a completely specular interface. 
In Figs. l4.5b )-f). the variation of the interface resistance as the function of vf2/vfi for 
different values of S is shown. In the absence of the potential barrier (vf2/vfi = 1), 
both Rq and Ri increase as the amount of diffuse scattering at the interface increase (S 
approaches 0), since only diffuse scattering is responsible for the resistance in this case. 
As the height of the barrier grows {vf2/vfi becoming > 1 or < 1), specular reflection 
starts to influence the interface resistance, and the values of both Rq and Ri are closer 
to each other for different values of S than they were in the absence of the potential 
barrier. Values of the potential barrier heights exist where the interface resistance is 
almost independent of the amount of diffuse scattering at the interface. These values 
are different for Rq and Ri. The percentage difference between Ri and Rq shows no 
variation with the barrier height in the case of completely diffuse interface (8 = 0), since 
VF2/VF1 only defines the scale of the interface resistance in this case. As the amount of 
diffuse scattering at the interface decreases, both — Rq)/Rq and (-Ri — Ro)/Ri show 
a strong dependence on the potential barrier height, reaching maximum values of several 
tens percent for small barriers. This fact is of particular interest, because small barrier 
heights are realized in Fe-Cr multilayers for minority electrons and Co-Cu multilayers 
for majority electrons. These systems will be discussed further below. 

4.1.3 Anticipation of the breakdown of the resistors-in-series model 

The results of the previous section show that the resistance of the interface deduced from 
the voltage drop measured far from the interface should not be interpreted as the 
interface resistance that would be obtained by measuring the voltage drop directly at the 
interface (Rq), as they may differ significantly. As seen from Figs. 14. lb ) and 14.31 the 
chemical potential varies within the electron mean-free path from the interface. Hence, 
in the multilayered system with layers thicknesses of the order or less than the mean-free 
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Figure 4.5: Summary of the results for the two-layered system: a) resistances Rq, Ri, 
and the difference Ri — Rq as the functions of the diffuse scattering S at the interface for 
vf2/vfi = 2, b) percentage differences Ri — Rq relative to _Ro and Ri as the functions of 
S for Vp2lvp\ = 2, c)-d) resistances Rq and -Ri as the functions of the electrons Fermi 
velocities ratio vf2/vfi for different values of the diffuse scattering S at the interface. 
In the plots a)-d), arbitrary units for the resistance are used. e)-f) percentage differences 
Ri — Rq relative to Rq and Ri as the functions of the electrons Fermi velocities for 
different values of S. 
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path, the chemical potential can be expected not to reach its asymptotic value within a 
layer. The breakdown of the resistors-in-series model can be anticipated, which means 
that while the resistance of the whole structure still can be found by adding the resis- 
tances due to the bulk of the layers and due to the interfaces, neither the resistance Ri 
nor Rq of the interface between the semi-infinite layers should be considered as the re- 
sistance of the interface between the layers which are thinner than the mean free path 
due to the scattering indigenous to them. In other words, even if one maintains the same 
interface between layers, its contribution to the resistance of a multilayered structures 
will depend on the thickness of the adjacent layers, or, more precisely, on the ratio of 
the thickness to the mean-free path. Its contribution will be constant, Ri, only if the 
interfaces are far enough from each other for the chemical potential to level off. 

In the following section, the prediction about the resistors-in-series model break- 
down and interface resistance dependence on the ratio of the layer thickness to the elec- 
tron mean-free path in this layer will be tested in the systems consisting of three and five 
layers. 

4.2 Three and five layers 

In this section, systems consisting of three and five layers of different metals are con- 
sidered. [56 1 Figure l4~6l shows the step-like potentials experienced by electrons in such 
systems. In the three-layered system ("Fig. 14.6b ')'). electrons move in a constant potential 
Vi = —mvpi/2 at z < —d and z > d, and V2 = —mvp2/'^ at —d < z < d. In the 
five-layered system (Fig. 14.6b )). electrons move in a constant potential Vi = —mvp-^^/2 
at z < —3d, —d < z < d, and z > 3d, and V2 = —mv'^p2/2 at —3d < z < —d and 
d < z < 3d. The structures with only two possible values of the electron Fermi velocity 
are chosen for simplicity; further below, when the realistic Fe-Cr and Co-Cu systems 
will be considered,this restriction will be removed, and the electron Fermi velocity will 
depend not only on the material but on the magnetization direction in the magnetic layer. 

Equations (IA.15I) - (IA.19I ) and (IA.20I) - (IA.28I) with the reflection coefficients (14.11) are 
solved in order to investigate how the chemical potential profile changes when the thick- 
ness of the layers changes and becomes less than the electron mean-free path in the 
layers. I find the resistance due to the interfaces as the function of the layers thickness 
relative to the electron mean-free path, for different values of diffuse scattering S at the 
interfaces, and investigate the deviation from the resistors-in-series model as the barrier 
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Figure 4.6: Step-like potential experienced by electrons in a) a three-layered system, b) 
a five-layered system. 



height Vp2/vpi changes, for different values of S. 

Figure 14771 shows the chemical potential profiles for three- and five-layered systems 
with S = 0.5; this represents interfaces where there is a 50% probability that the elec- 
trons are scattered diffusely and 50% specularly. In the case where the thicknesses 
of the inner layers are much larger than an electron mean free path in these layers 
(Fig. I4.7b ).b)'). the chemical potential approaches constant values in the inner layers. 
This means that the interfaces are independent of each other, and the interface resistance 
of the N-interface system is found to be N times the resistance of one single interface 
Ri, independently of the layer thicknesses as long as they are larger than the mean free 
path. Because of the diffuse scattering processes in the inner layers electrons lose their 
memory of the scattering at an interface and come to the next interface as if they were 
propagating in an infinite layer. In the second case, where the layer thicknesses are 
much smaller then the electron mean free path (Fig. l4.7b ').d')'). electrons retain a memory 
of the scattering from a previous interface; in this case the chemical potential does not 
approach a constant value, and different interfaces interact with each other, so that the 
interface resistance depends on the distance between the interfaces. 

Figure l4~8l summarizes the behavior of the interface resistance in the multilayers as 
the thickness of the layers, amount of diffuse scattering at the interfaces, and the height 
of the potential barriers experienced by electrons at the interfaces change. Arbitrary 
units for resistance are used, as in the Fig. 14.51 since the calculations are performed ir- 
respective of the particular metals. Figure R3k ) shows how the resistance due to two 
interfaces, R2, in the trilayer system depends on the thickness 2d of the middle layer 
relative to the electron mean free path I2 in this layer, for different amounts of diffuse 
scattering at the interfaces, and in the case where the electron Fermi velocity in the 
middle layer is twice smaller than the electron Fermi velocity in the outer layers, i.e.. 
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Figure 4.7: Chemical potential profiles for S = 0.5, normalized with respect to potential 
drop eEl within each layer: a) in the three-layered system with the middle layer thick- 
ness much larger than the electron mean-free path in the layer, b) in the five-layered 
system with the middle layers thicknesses much larger than the electron mean-free path 
in the layers, c) in the three-layered system with the middle layer thickness much smaller 
than the electron mean-free path in the layer, d) in the five-layered system with the mid- 
dle layers thicknesses much smaller than the electron mean-free path in the layers. 
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vpi = 2vf2- When the thickness of the middle layer is large (2d/l2 ^ 1), R2 is indepen- 
dent of 2d, and it is twice the resistance of one single interface Ri. As the ratio 2d/l2 
decreases, interface resistance in the trilayer system changes. For completely diffuse 
interfaces (S = 0) interface resistance decreases as 2d/l2 decreases. If a small amount 
of specular scattering is present at the interfaces, interface resistance at first increases 
significantly as 2d/l2 decreases, and then decreases abruptly when 2d/l2 becomes less 
than 0.2. As the amount of specular scattering at the interfaces grows (S reaches 1), the 
change in the interface resistance as 2d/l2 decreases becomes smaller. 

In the Fig. 14.8b ). the difference between the resistance due to the interfaces in the 
three-layered system R2 and twice the resistance of the single interface between two 
semi-infinite layers 2Ri relative to 2Ri as a function of the electron Fermi velocities 
ratio VF2/VF1 for different values of S at the interfaces and for 2d/l2 = 0.4 is shown. 
This difference represents the deviation from the resistors-in-series model. In the ab- 
sence of the barrier {vf2/vfi = 1), resistance exists due to diffuse scattering, and 
{R2 — 2Ri)/2Ri = for completely specular interfaces only. The deviation from 
the resistors-in-series model is the largest for the completely specular interfaces {S = 1) 
and VF2/VF1 = 1-4. The shape of the potential barrier (whether the electrons expe- 
rience a step up or step down at the interface between the first and the second layers) 
strongly influences the (i?2 — 2i?i)/2i?i difference, as the plots are asymmetrical around 
VF2/VF1 = 1, the ratio of the electron Fermi velocities that represents the absence of the 
barrier. 

Similarly, Fig. 14.85 ) shows the resistance due to four interfaces in a five-layered 
system as a function of the layers thickness 2d relative to the electron mean-free path 
in the second layer I2. Unlike in the three-layered system, the variation of the interface 
resistance with 2d/l2 is less pronounced for diffusive interfaces (S approaching 0), and 
more pronounced for specular interfaces {S approaching 1). Figure lOl i) shows the 
deviation from the resistors-in-series model in the five-layered system, the difference 
between R4 and 4_Ri relative to 4i?i for different S and 2d/l2 = 0.4. As the number 
of steps experienced by the electrons grows from 2 in the three-layered system to 4 in 
the five-layered system, the actual shape of the barrier seems to play less important role 
in defining the (Rn — NRi)/NRi difference as a function oi VF2/ vfi, the plots for 
(i?4 — 4i?i)/4i?i look symmetrical around vf2/vfi = 1- As in the three-layered system, 
the maximum deviation from the resistors-in-series model is obtained for completely 
specular interfaces; it reaches 15% ioxvF2/vFi = 1-6 and vf2/vfi = 0.6. 
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Figure 4.8: Summary of results for the three-layered and five-layered systems: a) resis- 
tance due to the interfaces R2 of the three-layered system as a function of the middle 
layer thickness to the electron mean-free path in this layer ratio 2d/l2 for different val- 
ues of diffuse scattering S at the interfaces and for the electron Fermi velocities ratio 
'VF2/VF1 = 0.5, b) percentage difference between R2 and twice the resistance of the 
single interface between two semi-infinite layers 2Ri relative to 2Ri as a function of 
the electron Fermi velocities ratio vf2/vfi for different values of S at the interfaces and 
for 2d/l2 = 0.4, c) resistance due to the interfaces R4, of the five-layered system as a 
function of 2d/l2 for different values of S at the interfaces and for vf2/vfi = 0.5, d) 
percentage difference between R4 and four times the resistance of the single interface 
between two semi-infinite layers 4i?i relative to ARi as a function of vf2/vfi for dif- 
ferent values of S at the interfaces and for 2d/l2 = 0.4, In the plots a) and c), arbitrary 
units for the resistance are used. 
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4.3 Co-Cu and Fe-Cr systems 



In this section the resistance due to the interfaces in the realistic two-, three-, and five- 
layered Fe-Cr and Co-Cu systems is studied. The following cases are considered: 1) a 
FMt-N single interface, 2) FMt-N-FMt(|)- a trilayer with parallel (anti-parallel) mag- 
netization directions in the ferromagnetic layers, and 3) FMt-N-FMt(|)-N-FMt - a five- 
layered system with parallel (anti-parallel) magnetization directions in the ferromagnetic 
layers (see Fig. 14.91) . In the case of the five-layered system the thicknesses of the inner 
layers are taken to be the same, and the magnetization of the middle ferromagnetic layer 
only is changed. Majority and minority electrons experience different potential steps at 
the interfaces between magnetic and non-magnetic metals (see Fig. l4.10b ')'). The ratio of 
the electron Fermi velocities in the Fe-Cr systems VF{Cr)/vF{Fe) is taken to be 0.837 
for majority electrons, and 1.003 for minority electrons. iBTll In the Co-Cu systems the 
ratio vp{Cu) /vf{Co) is calculated to be 1.05 for majority electrons, and 1.18 for mi- 
nority electrons (see appendix I A. 7. II) . Interface resistances are calculated by assuming 
majority and minority channels to be independent, [31 so that the total resistance is the 
parallel combination of the resistances of up electrons i?| and down electrons Ri : 

The step-like potentials experienced by electrons in the two-, three-, and five-layered 
Co-Cu and Fe-Cr structures are shown schematically in the Fig. 14.101 In order to find 
the resistances of up and down electrons and in each case a) through j) (see 
Fig. OUb . equations (im )- (in3l . (lAl3b - (lAl^ . and (lA^ - dA?^ were numerically 
solved with the reflection coefficients ( 14.11 ) and the appropriate parameters vpi, vf2, 
and vf3- I present the correspondence between R^, R^, type of the structure, equations 
that need to be solved, and parameters vpi, vf2, and vf3 entering these equations in 
the Appendix IA.7.21 In order to deduce the interface resistance from the jump of the 
chemical potential (Eq. ( 14. 31) ). the values of the electrons Fermi momentum in the non- 
magnetic layers are taken to be kpcu = l-36y4~^ for copper, and kFCr = S.8A^^ for 
chromium. [57| 

Figures l^.l H and 14. 121 summarize the behavior of the resistance due to the interfaces 
in the Co-Cu multilayers with different magnetization configurations as I change the 
amount of the diffuse scattering S at the interfaces and the thickness of the layers 2d 
relative to the electron mean-free path Imfp- Figure 14. lib ') shows the resistance _Ri 
measured far from the interface between two semi-infinite layers of Co and Cu, and that 
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Figure 4.9: a) FMt-N single interface, b) FMt-N-FMt(i) three-layered system, c) FMt- 
N-FMt(i)-N-FMt five-layered system. 
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Figure 4.10: Step-like potentials experienced by up (blue) and down (red) electrons 
in the a)-b) Co-Cu and Fe-Cr structures, c)-d) Co-Cu-Co and Fe-Cr-Fe structures with 
parallel magnetizations of the ferromagnetic layers, e)-f) Co-Cu-Co and Fe-Cr-Fe struc- 
tures with anti-parallel magnetizations of the ferromagnetic layers, g)-h) Co-Cu-Co-Cu- 
Co and Fe-Cr-Fe-Cr-Fe structures with parallel magnetizations, i)-j) Co-Cu-Co-Cu-Co 
and Fe-Cr-Fe-Cr-Fe structures with anti-parallel magnetizations. 
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measured right at the interface, Ro, as the functions of S. For diffuse interfaces (S is 
close to zero), Ri is larger than Rq. Both Ri and Rq decrease as the amount of diffuse 
scattering at the interface decreases (5* increases), but Ri decreases faster, and there 
exists a value of parameter S when Ri = Rq. The resistance Rq becomes larger than 
Ri for a specular interface. The relative difference between Ri and Rq (Fig. 14.1 lb )) is 
about 20% for completely diffuse interfaces; as the resistances decrease, this difference 
may reach more than 100%. 

Figures 14.1 lb )-f) show the resistances due to the interfaces in the multilayers R2 
for three-layered systems and R^ for five-layered systems as the functions of the layer 
thicknesses relative to the electron mean-free path for different amounts of diffuse scat- 
tering at the interfaces and different magnetization configurations (see Fig. 14.101) . In 
all four geometries, resistance due to the interfaces decreases as the amount of diffuse 
scattering decreases (S changes from to 1). Interface resistance is very small for com- 
pletely specular interfaces, due to the very small potential barrier height that the majority 
electrons experience at the interface between Co and Cu, so that the majority channel 
shunts the current. As seen from the Figs. 14.1 lb )-f). resistance due to the interfaces in 
the multilayers do not depend on the thickness of the layers 2d relative to the electron 
mean-free path Z^/p when 2d/lmfp is greater than 1, and they are found to be equal to 
the sum of resistances Ri of the corresponding number of the independent interfaces, 
2Ri or 4i?i. As the thickness of the layers becomes much smaller than the mean-free 
path, 2d/lmfp ^ 1, resistances R2 and R^ start to deviate from 2Ri or 4i?i, as predicted, 
indicating the breakdown of the resistors-in-series model. 

Figures I4.12b )-d) show the deviations of the actual resistances R2 and -R4 from the 
resistors-in-series model results, 2Ri or ARi, relative to 2Ri or 4i?i, as the functions of 
2d/lmfp for different amounts of diffuse scattering at the interfaces and different mag- 
netization configurations. These deviations are not monotonic with 2d/lmfp- The differ- 
ence between R^ and NRi is especially large for extreme interfaces, almost completely 
specular, S close to 1, or completely diffuse, S close to 0. (In the three-layered systems 
the plot for S" = 1 is not shown, as it has a much larger scale than other plots.) It reaches 
10-15% in the three-layered systems, 20% in the five-layered systems with the paral- 
lel magnetization directions, and 30% in the five-layered systems with the anti-parallel 
magnetization directions for completely specular interfaces. 

Finally, Figs. I4.1^ )-f) show the difference between actual resistance due to inter- 
faces Rm (N =1,2, and 4) and the resistance that would be obtained assuming interfaces 
to be independent of each other as a function of the amount of diffuse scattering S at 
the interfaces. The thickness 2d of the non-magnetic layers is taken to be small com- 
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pared to the electron mean free path in the layers {2d/lmfp=0A). The actual resistance 
is compared with both the resistance of an independent interface measured right at the 
interface, NRq, and that measured far from the interface, NRi. In all geometries, both 
{Rn — NRi)/NRi and (Rn — A^-Ro) /NRq are significant for either completely diffuse 
or completely specular interfaces. For the realistic interfaces where an equal amount 
of specular and diffuse scattering is present, S = 0.5, the difference between the inter- 
face resistance one could measure and the resistance of independent interfaces measured 
right at the interfaces, NRq, reaches 10%, which is comparable with the experimental 
error for the interface resistance. Il23ll 

Figures 14.131 and 14.141 show the same information as the figures 14.1 II and 14.121 for 
Fe-Cr systems. Similarly to the resistances in the Co-Cu systems, Ri is larger than Rq at 
small S, both Rq and Ri decrease as the interface becomes more specular (S approach- 
ing 1), Ri decreasing faster (see Fig. I4.13h ')'). The difference with Co-Cu systems is 
that Ri becomes equal to Rq for S closer to 1. The relative difference between _Ro and 
Ri is between 12 and 18% for S up to 0.9.(Fig. 14.13b ')'). after that the difference grows 
significantly (not shown). Resistances R2 and R4 for different magnetization configura- 
tions decrease as the amount of diffuse scattering at the interfaces decreases. In the case 
of the parallel magnetization directions both in the three- and five-layered structures, 
resistances of the completely specular interfaces are not shown, because they are small 
compared to the resistances for S 1, due to the very small potential barrier experi- 
enced by minority electrons at the interface between Fe and Cr (smaller than that for 
majority electrons at the interface between Co and Cu). Resistance due to the interfaces 
depend on the thickness of the layers relative to the electron mean-free path, 2d/lmfp, 
as it becomes of the order or less then 1 (Fig. I4.13b ')-f)'). As it is for the Co-Cu sys- 
tems, the deviation from the resistors-in-series model for Fe-Cr structures is substantial 
for extreme interfaces, with S close to or 1 (Fig. l4.l4b )-d)'). but for the interfaces with 
5 = 0.5, the difference between R^ and NRq reaches 20% for the systems with parallel 
magnetic configuration and more then 15% for the systems with anti-parallel magnetic 
configuration (Fig. l4.13'fe ')'). This exceeds experimental error, which is about 10% for the 
interface resistance and about 15% for the bulk resistance ll23l . 



60 



Co-Cu 




AR2 ,fQ.n? 



parallel 



c) 



2 

1.5 

1 : 

0.5 




s=o 

S=0.1 

S=0.2 
S=0.3 
S=0.4 
S=0.5 
S=0.6 
S=0.7 
S=0.8 
S=0.9 



42dAmfp5 



''2dAmfp^ 



AR2 ,fD.n? 



anti -parallel 



d) 



1.5 
1 



S=0 

S=0.1 
S=0.2 



0.5 — 

— 



S=0 
S=0 
S=0 
S=0 
S=0 
S=l 



2d A, 



■mfp 




AR4 ,fnir? 



anti -parallel 



f) 



s=o 

S=0.1 
S=0 .2 



S=0.3 
S=0.4 
S=0.5 
S=0.6 
S=0.7 
S=0.8 
S=0.9 
S=l 



'2dAmfp5 



Figure 4.11: Resistances due to the interfaces in the Co-Cu systems: a) sheet resistance 
measured far from the interface, ARi, and that measured at the interface, ARq, in the 
two-layered systems as the functions of the amount of diffuse scattering S at the inter- 
face, b) the difference between Ri and Rq relative to Ri in the two-layered systems as 
a function of S, c)-d) resistance R2 as a function of the inner layer thickness relative to 
the mean-free path in the three-layered systems for different values of the diffuse scat- 
tering S and different magnetization configurations, e)-f) resistance R4 as a function of 
the inner layers thicknesses relative to the mean-free path in the five-layered systems for 
different values of S and different magnetization configurations. 
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Figure 4.12: Deviations from the resistors-in-series model in the Co-Cu systems: a)-b) 
the difference between the resistance R2 and twice the resistance of a single interface 
2Ri relative to 2Ri as a function of the inner layer thickness 2d relative to the mean-free 
path Imfp in the three-layered system for different values of S and different magnetiza- 
tion configurations, c)-d) the difference between the resistance i?4 and four times the 
resistance of a single interface, ARi, relative to ARi as a function of 2d relative to Imfp 
in the five-layered system for different values of S and different magnetization configu- 
rations, e) the difference between the resistances Rf^, with = 1, 2, 4 and times the 
resistance i?o, NRq, of a single independent interface measured right at the interface, 
relative to NRq, in the multilayers with 2d -C Imfp as a function of S for different mag- 
netization configurations, f) the difference between the resistances R^, with = 2, 4 
and times the resistance Ri of a single independent interface measured far from the 
interface, NRi, relative to NRi, in the multilayers with 2d <ti Imfp as a function of S 
for different magnetization configurations. 
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Figure 4.13: Resistances due to the interfaces in the Fe-Cr systems: a) sheet resistance 
measured far from the interface, AR\, and that measured at the interface, Ai?o, in the 
two-layered systems as the functions of the amount of diffuse scattering S at the in- 
terface, b) the difference between R\ and i?o relative to i?o in the two-layered systems 
as a function of 5" (note the different S-axis limit), c)-d) resistance i?2 as a function of 
the inner layer thickness relative to the mean-free path in the three-layered systems for 
different values of the diffuse scattering S and different magnetization configurations, 
e)-f) resistance i?4 as a function of the inner layers thicknesses relative to the mean-free 
path in the five-layered systems for different values of S and different magnetization 
configurations. 
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Figure 4.14: Deviations from the resistors-in-series model in the Fe-Cr systems: a)-b) 
the difference between the resistance R2 and twice the resistance of a single interface 
2Ri relative to 2Ri as a function of the inner layer thickness 2d relative to the mean-free 
path Imfp in the three-layered system for different values of S and different magnetiza- 
tion configurations, c)-d) the difference between the resistance i?4 and four times the 
resistance of a single interface, 4i?i, relative to 4i?i as a function of 2d relative to Imfp 
in the five-layered system for different values of 5* and different magnetization config- 
urations, e) the difference between the resistances R^, with = 1,2, A and N times 
the resistance Rq of a single independent interface measured right at the interface, NRq, 
relative to NRq, in the multilayers with 2d <^ l^fp as a function of S for different mag- 
netization configurations, f) the difference between the resistances R^, with N = 2, A 
and N times the resistance Ri of a single independent interface measured far from the 
interface, NRi, relative to NRi, in the multilayers with 2d -C Imfp as a function of S 
for different magnetization configurations. 
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Chapter 5 

Spin-torques - theory 



In the following two chapters, a mechanism of the magnetization switching of magnetic 
multilayers driven by a spin-polarized current is studied. First, I present a spin trans- 
fer model in which the equation of motion of the spin accumulation is solved in order 
to derive the torque acting on the background magnetization. Next, a system of two 
thick magnetic layers separated by a non-magnetic spacer is considered. The analytical 
expressions for the spin-current and spin-accumulation distributions in the system, and 
spin-torque and effective field acting on the local magnetization in the ferromagnetic 
layers are obtained, and I discuss the significance of the result for the spin-current at 
the interface between layers in the system with magnetization directions close to anti- 
parallel. Finally, the correction to the CPP magnetoresistance due to the interaction of 
the electron spins and background magnetization is calculated. In the next chapter, I 
present the results for a realistic multilayered structure. 



5.1 Review of formalism 

The linear response of the current to the electric field in the x-direction for diffusive 
transport, Eq. (12.261) . can be written in a spinor form as 1 361 

j(x) = CE{x)-D—, (5.1) 

where E{x) is the electric field, j, C, D, and h are the 2x2 matrices representing the 
current, the conductivity, the diffusion constant, and the charge and spin accumulations 
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at a given position. The diffusion constant and the conductivity are related via the Ein- 
stein relation C = e^N{eF)D for a degenerate metal, where N{eF) is the density of 
states at the Fermi level. In general, these matrices can be expressed in terms of the 
Pauli spin matrix cr as 

' C = CoI + (T-C, 
D = DqI + a - B, 

N = Noi + fT-'N, (5.2) 
n = UqI + cr • m, 
j = jj + cr ■ j^, 

where 2no is the charge accumulation, 2no = + n^, and m is the spin-accumulation, 
2|m| = n| — nj, n| is the number of electrons with the spins parallel to the princi- 
pal quantization axis, is the number of electrons with the spins anti-parallel to the 
principal quantization axis, and, similarly, electric current 2je = Jt + ii' magnetization 
current 2|jm| = Jt — Ji- Inserting Eqs. (15.21) in Eq. (15.11) . the electric current je and 
magnetization current can be written as 

je = Re(Trj) = 2CoEix) - 2Do^ " 2° " (5-3) 

= ReTr(fTj) = 2CE{x) - 2D^ - 2/^o^, (5.4) 

where the units e = fiB = I are chosen for the notation convenience. The first term 
on the right hand side of Eq. (15.41 ) is the contribution to the spin polarized current from 
the electric field, the third term is the contribution from the spin accumulation, which is 
present in a magnetically inhomogeneous structure. The second term, the contribution 
coming from the charge accumulation, will be neglected, since its characteristic length- 
scale, of the order of several angstroms, is of the order of magnitude smaller than other 
relevant lengthscales. As will be shown below, in magnetic multilayers the contribution 
of the third term to the spin current can dominate over the first. Comparing Eq. (15.31 ) 
with the expression (12.261) for the electric current, the parameters Cq and Dq can be 
identified with Cq = ^cr = Dq = \D = \vIt = IvpXmfp, where A„/p = vpT 

is the electron mean-free path. 

For a transition metal ferromagnet, the spin-polarization parameter for the conduc- 
tivity P, the diffusion constant /?', and for the density of states P" can be defined as 
C = pCoMd, D = P'DoMd, and N = P"NoMd, where is the unit vector to 
represent the direction of the local (background) magnetization. Similarly to the spin- 
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polarization parameters for the interface, 7, 7', and 7", introduced in the Appendix lB.il 
spin-polarization parameters for the bulk are related as = {[3 — [3") / {1 — (3 (3") (com- 
pare with Eq. (IB ■341) ). The experimental value for [3 can be found in the Ref. ll23ll . and 
the value of [3" can be found from the bulk densities of states of up and down electrons, 
and N^, as [3" = {N^ - Ni)/{N^ + N^). Inserting the expressions for C and D in 
Eqs. (15.31) and (15.41) . and eliminating the electric field and charge density, the following 
expression for the spin-current is obtained: 



dm 

dx 



- /3/?'Md(M, 



dm. 

dx 



(5.5) 



The equation of motion for the spin accumulation may be obtained by generalizing 
the continuity equation (12.241) to the systems with both spin-up and spin-down electrons, 
and including the possibility of spin-flip scattering. The continuity equation for the spin- 
up electrons density (spin-down electrons density n^) has the form 



dt dx 



where r^/ is the spin-flip relaxation time of the conduction electrons, and, for the spin- 
accumulation, 

d\m\ ^ d\]m\ _ 2|m| 



dt dx r. 



Introducing the term J/hm x (see Chap.[T]), the equation of motion for m may be 
written as 

dm dim J ,^ 2m 

^ + 1^ + T^n X Mrf = . (5.6) 

dt dx n Tgf 

Upon placing the expression for the magnetization current Eq. (15.51) in Eq. (15.61) one 
finds the following equation of motion for the spin-accumulation vector: 

1 dm d^m , ^ (9^m\ m m x 



2Do dt dx^ ' ^'^^ V'^" dx^ / A2 ' A2 



^ (/3jeMrf), (5.7) 



2Dq dx 

where A^/ = J D^Tsf and Aj = J^UDs^j J . The term on the right hand side of the 
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time dependent diffusion equation Eq. ( 15 .71 ) for the spin accumulation is the source 
term; it is this term that drives the accumulation. |j59| Following the conventional treat- 
ment of magnetic multilayers for current perpendicular to the plane of the layers, one 
can assume that the magnetization is uniform throughout a layer, and changes discon- 
tinuously between the layers. 11601 i6 1 If one looks for the steady-state solutions of the 
equation (15.71) . the first term on the left-hand side is zero, and the electric current je is 
constant throughout the multilayer; therefore one can see that the source term is confined 
to interfaces between the layers, and can be taken into account by appropriate bound- 
ary conditions. The boundary conditions at the diffusive interfaces between two layers 
with noncoUinear magnetizations are addressed in detail in Appendix lB.il It is shown 
that while the spin- accumulation experiences a jump at the interfaces, proportional to 
the interface resistance, the spin-current is continuous across the diffusive interface, and 
no spin-torque is created at the interface. The specular scattering at the interfaces can 
also be taken into account within the semiclassical framework presented in this work, 
as it is done, for example, in Chap. |3 this involves considering each k vector sepa- 
rately which for a noncoUinear multilayer is a cumbersome problem, and it will not be 
addressed here. Below, the equations for the spin-accumulation and the effect of the 
spin- accumulation on the local magnetization in the bulk of the layers are discussed. 

The equation of motion for the local magnetization is the Landau-Lifshitz- 
Gilbert equation 

^ = -7oM, X (H, + Jm) + aM, x (5.8) 

at at 

where 70 is the gyromagnetic ratio, Hg is the magnetic field including the contributions 
from the external field, anisotropy and magnetostatic field, the additional effective field 
Jm is due to coupling between the local moments and the spin accumulation, and the 
last term is the Gilbert damping term. Note that the characteristic time scale for the local 
moment is 'jQ^He, of the order of nanoseconds for a magnetic field of 0.1 T, while the 
time scale of the spin accumulation is of the order of r^j and h/J, i.e., of the order of 
picoseconds. Therefore, the background magnetization can be assumed to be fixed on 
the time scale of the spin accumulation, and the steady state solution of the equation (15.71) 
can be considered. 

By separating the spin accumulation into longitudinal (parallel to the local moment) 
and transverse (perpendicular to the local moment) modes, m|| and in±, and looking for 
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the steady state solutions, equation (15.71 ) can be written as 



(9^mii m 



where Kdi = a/I - I3l3'\sf, and 

9^m_|_ m_|_ m_|_ x 



0. (5.10) 



Note that the terms longitudinal and transverse refer to the magnetization in the indi- 
vidual layers, i.e., they are locally defined and have no global meaning throughout a 
noncoUinear multilayer. 

The longitudinal accumulation m|| decays at the length scale of the spin diffusion 
length Xsdi, while the transverse spin accumulation in± decays as A j if one assumes 
Aj -C As/. This assumption is valid for cobalt, for example, where the spin-diffusion 
length Xsdi has been measured to be about 60 nm, /3 is 0.5, A^/p is about 6 nm, lE^ 
vp = 0.6 ■ 10^ m/s, J = 0.3 eV.(W\ I estimate jS' to be about 0.9 using the bulk densi- 
ties of states for up and down electrons of 2.42 states/(atom Ry) and 13.42 states/(atom 
Ry) ll6Tl . so that Xgf = 80 nm, Aj = 2 nm, Xj -C As/, and the transverse spin accumu- 
lation has a much shorter length scale compared to the longitudinal one. However, for 
permalloy, As^; has been measured to be about 5 nm, f3 is 0.7, lE^ (3' is estimated to be 
about 0.95 using the bulk densities of states for up and down electrons of 2.4 states/Ry 
and 15 states/Ry ll^l^ . so that As/ is about 9 nm. Taking the typical diffusion constant 
of a metal to be Dq = 10~^m^/s, and J = 0.1 eV, [41 1, Xj can be estimated to be 4 nm, 
comparable with As/. The majority of the results presented in this work are obtained in 
the limit Xj <^ As/, so that they are applicable to cobalt. The structure consisting of two 
semi-infinite permalloy layers divided by a nonmagnetic spacer is considered separately 
below. 

As seen from Eq. (15.81 ). the longitudinal spin accumulation has no effect on the local 
moment; therefore, one can re-write Eq. (15.81) in terms of the transverse spin accumula- 
tion only by replacing m by : 

^ = -7oM, X (H, + Jm^) + aM, x (5.11) 
at at 

In order to describe the transverse accumulation, it is convenient to introduce an aux- 
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iliary vector A such that Jm± = A x M^. If one considers a system with two 
noncollinear ferromagnetic layers, the spin accumulation in one layer depends on the 
orientation of the other. Supposing that the above equation is used for the layer Fl (see 
Fig. ll.lll) . the local magnetization of this layer is labeled as M^^'*, and the magnetiza- 
tion of the other layer is labeled as M^^\ Without loss of generality, the vector A can 

(2) 

be written as a linear combination of the vector , and the vector perpendicular to 

(2) (2) (2) (1) 

M); ' , for example, A = aM)^ ' — 6M)^ ' x M)^ ' , so that the two components of the 
accumulation in the plane perpendicular to M^^^ are written as 

Jm± = aMj^ X M^^^ + bM^J'^ x (mJ^ x mJ^), (5.12) 

where a and b are determined by geometric details of the multilayer. By placing this 
form of the accumulation in the equation of motion for the background magnetization, 
Eq. (15.111) . one finds the following equation of motion for M^j^-*: 

^ = -7oM(^) X (H, + 6M?))-7oaMi^) x (M?) x M«) +aM, x (5.13) 



Thus the transverse spin accumulation produces two effects simultaneously: the term 
bM.'^l^ X M^^^ is the torque due to an "effective field" bMf\ and the other is aM^^'' x 
(Mf^ X M^^^) is called the "spin torque" predicted in the Refs. 10 El S COl Hi | . The 
first term produces a precessional motion about M^j"^^; in this sense it acts as if the spin 
current creates a magnetic field on M^j"'^'' (see Fig. 15.11) . The second term acts so as 
to increase or decrease the angle between M.^J^^ and M^f'*; also, it acts so as to assist 
or oppose the damping term in Eq. ( 15.81) . Note that both the spin-torque and effective 
field terms appear on an equal footing in the equation ( I5.13I ). as both are related to the 
transverse spin accumulation. 

The torque transmitted by the current to the background may also be determined 
from the following considerations. The torque r acting on the background angular mo- 
mentum is the time derivative of M^, r = dM.d / dt. The total angular momentum of 
the system, which is the sum of the electron angular momentum m and the background 
angular momentum is conserved, so that m + M^^ = const. The torque acting on 
the background magnetization can be written in terms of the electron spin accumulation 
and spin current as 

dm dm djm 

dt dt dx ' 
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Figure 5.1: Direction of the effective field (red vector) and spin-torque (blue vector) 
acting on a thin FM layer background magnetization M^^"* = e^, if the magnetization of 
the thick layer is M.)]^ = cos 6ez — sin 6ey. 

so that the torque transmitted by a steady-state current is given by the gradient of the spin 
current. Integrating this over a layer with the thickness d which absorbs the momentum, 
one finds 



The total torque absorbed by the layer is the difference between the spin current values 
at the boundaries of the layer. 

5.2 Two ferromagnetic layers 

In order to illustrate the method developed in Sec. l5.1l to determine the effective field and 
spin-torque acting on the ferromagnetic layer, I choose a simplified system where I can 
analytically derive the spin accumulation, spin-current, effective field, and spin-torque. 
I consider a system consisting of two thick ferromagnetic layers with the magnetizations 
M^^'' = cos + sin 6ey and M^^'* = cos 6*62 — sin 6ey, where x is the direction of the 
electric current (see Fig. l5.2l) . separated by a nonmagnetic spacer layer so that there is 
no exchange coupling between ferromagnetic layers, and backed by nonmagnetic layers, 
or leads. The nonmagnetic spacer layer thickness is considered to be small compared 
to the spin-diffusion length in the nonmagnetic metal (Xf^i ~ 600 nm), so that both 
spin accumulation and spin current are constant in the spacer, and I will use both terms 
"spacer" and "interface between magnetic layers" interchangeably. The thickness of 




(5.14) 
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Figure 5.2: System of two thick ferromagnetic layers separated by a spacer layer 
and backed by nonmagnetic layers, or leads. Magnetization of the first (right) layer 
is M^^'' = 0086*62 + sin6'ey, magnetization of the left (second) layer is M^^"* = 
cos 6ez — sin 9ey. (x, y, z) is the global coordinate system, (a;, |/, z) and (x, y, z) are the 
local coordinate systems connected to the local magnetizations in the first and second 
layers. 



the ferromagnetic layers is large enough so that one can neglect the reflections from 
the outer boundaries of the layers. The difference between this system and the system 
considered in the Ref. I36J is that in the latter, the second layer is assumed to be half 
metallic, so that the current is fully spin-polarized, while here this approximation is 
abandoned. 

I solve Eqs. ( 15.91) and ( 15.101) to find the longitudinal and transverse spin accumulation 
in both ferromagnetic layers, and use Eq. ( 15.51) to find the spin-current. In the local 
coordinate system (x, y, z) connected to the local magnetization of the right layer, the 
spin- accumulation and spin current take the following form: 




(5.15) 
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Jm,z 



Jm,x 



-^sdl 



2Do(l-/3/30 



where Gi and G2 are constants of integration, and 



1-i 
V2X'j 



(5.16) 



(5.17) 



when Aj <^ A^/. Similar expressions for the local spin- accumulation and spin-current 
are written for the second layer. 



To determine the constants of integration, the boundary conditions at the interface 
between two ferromagnetic layers have to be invoked. The boundary conditions are dis- 
cussed in detail in Appendix lB.il At a diffusive interface, the spin-current is conserved, 
and the spin-accumulations to the left and to the right of the interface in the global 
(x, y, z) coordinate system are related as follows: 

m,(0-) - m,(0+) = 2j^,.(0)r^ 
m,(O-)-m,(O+) = 2j„,,(O)(cos2 0+ ^ 



(1-77') y V2D0 



[ "^.(0-) - m,(0+) = -7JeCOsgr ^^^_^^^,^^^ + 2j„,,(0) (sin^ ^ + 



(1-77')/ v^Do' 



where r is proportional to the interface resistance ARj (see Eq. ( IB.35I) ). 7 and 7' are the 
spin polarization parameters for the conductivity and the diffusion constant at the inter- 
face (see Eq. (IB. 51) ). Note that the spin-accumulation is conserved across the interface if 
there is no diffuse scattering at the interface {ARi = 0). Vectors in the local and global 
coordinate systems are related via the rotation matrix (see Eq. (IB.20l) '). and the constants 
of integration are 

^ ^ X -^'^i^ + ^) + '-^^o^'d 

' V2XjJefiB cos2^(l + Y^) + Asin2 0(l + f) 
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^ ^ Pje hal (l+ftA)sin2g _ 

' 2v/2AjJe/iB cos2^(l + ^) + Asin2^(l+f)^ ^' ^ 

where A = (1 - /?/5')Aj/v^A,rfz, t = ^/f3{l - 77'), ^ = V2ARje'N^{eF)Do{l - 
7"^)/Aj(l — 77'), ARi is the resistance of the interface between two ferromagnetic 
layers, /?, P' are spin-polarization parameters for conductivity and diffusion constant 
in the bulk of the layers, 7, 7', 7" are spin-polarization parameters for conductivity, 
diffusion constant, and electron density of states at the interface; Nq (ep) is the density 
of states at the interface at Fermi energy (see Appendix IB .11) . The expressions for the 
spin- accumulation in the right layer are 



m 

m. 
m 



«=2e-v^^^eG,[cos^ + sin^] 

« = _2e -^^jReG2[cos ^ - sin ^] cos^ + de"^ sin^ (5.20) 
41) = 2e"^ReG2[cos ^ - sin ^] sin^ + de"^ cos^. 

In the left layer they are 

7m^\x) = —m^^\—x). 

The spin-current in the layers can be found from Eq. (15.51) . The expressions for the 
components of are complicated, and I will only present the values of the x—, y—, 
and 2;— components of at x = 0, in the nonmagnetic spacer, in the absence of the 
scattering at the interface. Due to the symmetry of the problem, only the 2-component 
is non zero, so that 

cos^ y + A sni 6 

where A = a/1 — PP'Xj/\/2\sf, and it is of the order of 0.02 for cobalt. The spin 
current at the interface reaches its maximum value of 

3m,max[x = 0) = ^ (5.22) 
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when the angle between the local magnetizations is 26*, where 

1 _ Q \ 

cos2r = -, (5.23) 

1 — A 

and 29* is close to vr when A ^ 1. The magnitude of the spin-current is enhanced by a 
large factor of A^^ compared to the bare spin-current /3je cos 9. This enhancement comes 
from the interplay between longitudinal and transverse accumulations; it is the result of 
the global nature of the spin-current even though the transverse components of the spin 
current and accumulation are absorbed within a region of several Aj of the interface. 
While the total spin torque acting on a symmetric two-layered structure is zero, the 
result of Eq. (15.211) can serve as an indication that in the non-symmetric structure when 
one layer is pinned and the other free; the spin-torque acting on the free layer, and, 
hence, the angular momentum transferred to it, exceeds the transverse component of the 
bare portion of the incoming spin-current by 1/2 A. This prediction will be tested in the 
next chapter, where the non-symmetric three-layered system is discussed. 

The spin-torque oi and the effective field bi acting on the ferromagnetic layer can 
be found from the following considerations. From Eq. (15.151) . the transverse spin- 
accumulation in the right layer can we written as 

mf^ = m^^^ ex+m.il'^ By = 2Re{G2exp{—x/l+))e^+2lm{G2exp{—x/l+))ey, (5.24) 

where e^. and are unit vectors in the local coordinate system. By noticing that G2 can 
be written as G2 sin 29, M^^^ x M^^^ = - sin 29es, M^^^ x (M^^ x M^^^) = - sin 29ey, 
one can write m.^^^ as 

mi^^ = -2Re((5'2exp(-a;//+))My^ X M^^^ 

- 2lm{G2exp{-x/l+))M'-^'^ x (mJ^ x mJ^). (5.25) 

One can see that the form of the transverse spin-accumulation given by Eq. (15.251) is 
exactly the form used in the definition of a and b (Eq. (l5.12l) '): 

Jm± = aM^^ X M^^^ + bM^j^ x (mJ^ x M^^^). (5.26) 

To obtain the coefficients a and b, one averages the spin-accumulation (15.251) over the 
thickness ti;' of the right layer. The spin-torque and effective field per unit length acting 



75 



on the ferromagnetic layer take the form 



nal 1 - cos(^)e v^^. (1 + ftX) 



ai = ^ X —pi X 



2V2\jei^B tF/iV2\j) cos2^(l + ^) + Asin2^(l ' ' 



(5.27) 



' 2V2XjefiB W(v^Aj) cos2^^(l + ^)+Asin2^(l + f)' ^ ■ 

Note that in order to get the true torque and field acting on the ferromagnetic layers, the 
coefficients a and b has to be multiplied by sin 29 (see Eq. (I5.26I) ). 

Both spin-torque a sin 26* and effective field 6sin2^^ reach their maximum value 
when 9 = 9** where 

cos2r* = -^^ (5.29) 
1 + A ^ ^ 

in the absence of the interface resistance; we see that 9** is close to 90°. 

In the second layer, the spin-torque 02 is equal to — ai, and the effective field 62 is 
equal to —61 due to the symmetry of the problem. 

In Fig. 15.31 1 present the spin accumulation and spin current distribution in the FM- 
Sp-FM structure with Aj = 4 nm, Xsdi = 60 nm in the presence of the interface resis- 
tance ARj = 0.5 ffi-m^; the diffusion constant is taken to be 10^^ m^/s. Both m 
and jm are constant in the spacer, so I don't show the spacer in these plots. As follows 
from Eq. (15.201) . the characteristic length scale of the x-component of the accumula- 
tion, perpendicular to the plane of the layer magnetizations, is Aj. Global y- and z- 
components of spin accumulation are defined by both transverse and longitudinal local 
components of the accumulation, so for them, two length scales may be distinguished: 
one of the order of A j for a more rapid change of the spin accumulation, and the longer 
one, of the order of the spin-diffusion length Xsdi- Only the z-component of the spin 
accumulation is discontinuous in the presence of the interface resistance due to the sym- 
metry of the problem. Far from the interfaces, x ^ Xj and x <^ —Xj, the spin current 
is coUinear with background magnetization M^^-* or M^^-*, approaching its bare values 
jm = Pje^li^ or jm = /3jeM^^\ As foUows from the equation (I5.21I ). a large ampli- 
fication of z-component of the spin-current occurs within a distance of several Aj from 
the interface. 



In Fig. 15.41 1 show the spin-torque and effective field per unit length acting on the 
first ferromagnetic layer as a function of the thickness of the layers tp relative to V2Xj. 
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Figure 5.3: x-, y-, and 2;-components of the spin- accumulation 
T^/i/3je/V^^.jJ){fi'0'o/efJ,B) (a-c) and spin-current jm//?je (d-f) distribution in 
FM-Sp-FM system in the presence of interface resistance for A j = 4 nm, Xg^i = 60 nm 
and different angles 9. 
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Figure 5.4: Spin-torque a/ Pje{hal/eiXB) and (b) asm29 / Pje{hal/ efis), (c) effective 
field b//3je{hal/efiB) and (d) b sin 26 //3je{hal/ 6^13) acting on the FM layer in FM-Sp- 
FM system as a function of tp/\^Xj in the presence of interface resistance for Aj = 4 
nm, Xsdi = 60 nm and different angles 9. 



One can see that while the torque and field terms a and b (Fig. 15 .4^ ) and c)) are largest 
for 9 = 90°, they don't act on the background magnetization, since sin 29 = 0. The true 
torque and field a sin 29 and b sin 29 (Fig. l5.4b ) and d)) are largest for 9 close to 90° (see 
Eq. (15.291) ). It is worth noting that while the effective field decreases much faster then 
the spin torque as the thickness of the layers increases, at their maximum, they have the 
same magnitude, so they indeed have to appear on an equal footing in the equation of 
motion for the local magnetization. 
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5.3 Correction to CPP resistance 



The normalized angular dependence of the resistance can be defined as ll43l 

_ Ri29) - R{0) 
" ~ R{27r) - R{Oy ^ ^ ^ 

and the data was fit to 

_ l-cos^(g) 
1 + xcos^(t^) 

I calculate the angular dependence of the CPP resistance based on the spin currents 
that are found using the diffusion equation for the spin accumulation in noncoUinear 
structures as follows. 

The resistivity of a system is a proportionality coefficient between electric field and 
electric current, E = pj^,. From Eq. (I5.3I ). electric current can be written in terms of the 
longitudinal spin accumulation m|| as 

je = 2CoE-2Do(3'^, (5.32) 
ox 

where in the right layer m|| = = Gi exp (—x/\sdi), in the left layer m|| = mj = 
— Giexp {x/Xsdi) (see Eq. (15.151) ). Gi is defined by Eq. (15.181) . The resistivity of the 
first (right) and the second (left) layers then takes the form: 

Pl,2 - ^^jr • (3.ii) 

The sheet resistance of the whole system can be found by integrating the resistivity over 
the x-coordinate from —tp to tp, and adding the interface resistance ARj; I find 



AR(f)) = ^ + /?/3-A (1 + f) (1 - 770 sin^ 9 + f{l - 7//?) cos^ 9 \ 
Co V (1-/3/3') (1-77' + rA)cos2^ + A(l + f) (1-77') sin^^j 

(5.34) 

where 9 is half of the angle between magnetizations. The parameter x then takes the 
following form: 

1 f 

X = -1 + , + , -w. (5-35) 

A(l + r) (l+r)(l-77') 
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where A = (1 - (5P')\j/{^Kai), f = V2ARie''Nl{eF)D,{l - - 77'), e 

is the electron charge, Nq is the density of states at the interface, Dq is the diffusion 
constant in the bulk, 7, 7', 7" are the spin-polarization parameters for the conductivity, 
diffusion constant, and density of states at the interface (see Appendix lB.il) . For cobalt, 
I estimate x to be about 60 if I take Aj = 2 nm, Xgdi = 60 nm, Dq = 10^'^m^/s, P = 0.5, 
7 = 0.75, = 2.42 states/(atom Ry), A^^ = 13.42 states/(atom Ry) for the bulk, 
iV| =1.83 states/(atom Ry), = 17.75 states/(atom Ry) for the interface (with Cu), 
so that /?' = 0.9, 7' = 0.97, A^o(ei7') = 6.5-102^ states/m^eV. To my knowledge, there 
is no experimental data for the value of x for cobalt, so my result can not be compared 
with the experimental one. Below, I present a calculation for permalloy, for which an 
experimental value of x exists. 1 43J 



5.4 Permalloy 

For permalloy, the approximation Aj ^ A^/ is not valid, so one has to use the exact 
expression for /+ (see Eq. (15.171) ). 



1 



1 i c — id 

(5.36) 



1+ \ ^If v^Aj 



where c = ^{Xj/Kf^ + ^(Xj/Kf)^ + 1, d = ^-{Xj/KfY + ^{Xj/Kf)^ + 1. 

Note that for cobalt c ~ ~ 1. In the absence of the interface resistance, I obtain 
the following analytical expressions for spin-accumulation distribution, spin-torque, and 
effective field: 

^{1) = 2e"^ [ReGa cos ^ - ImGz sin ^] 

m« = -2e~^ [ReGs cos ^ + ImGa sin ^] cos 6 + de"^ sin 6 
= 2e"^[ReG2Cos^ + ImG'2sin ^] sin^ + de"^ cos^, 



J l-cos(J^)e"^^ 



C2 + rf2 tF/{V2X 
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J . . sin(^^)e 

-(cReGs-rflmGa) ^ ' 



J 



J ^ _ 1 - cos(^^)e 

-XdReG2 + c\mG2) ^ , (5.38) 



c2 + ci2^ tp/iV^X 



J 



where 



^ XjJ efXB ic^ + d^)cos^e + 2cXsm^e ^ ' ' 



V2XjJefiB (c2 + rf2) cos2 ^ + 2cA sin^ ^ ^ ^ 

The parameter x (Eq. l5.3ll) for permalloy takes the form 

c^ + ci^ 

X^-l.^, (5.4,) 



where 



]l{Xj/x,fy + ^{Xj/Xsf)^ + i 



+ ^2 ^ 2v^(Aj/A,j)4 + 1 = 2^{l-(3mXj/XsdiY + l, 

^_ (1-/3/30 Aj 
-\/2 Asrf/ 

Taking Aj to be 4 nm, As^^ to be 5 nm, /3 = 0.7, /?' = 0.95 (see Sec. 15. IL I estimate 
X to be about 3.8. This number has to be compared with the experimentally obtained 
X = 1.17. [43 1 While the theoretical value of x is more than three times larger than 
the experimental one, it is calculated without taking into account the resistance of the 
interface between permalloy layers; the presence of he interface resistance leads to a 
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smaller value of x (see Eq. ( I5.35I) ). 



Chapter 6 

Spin-torques in the thin FM - thick FM 
- NM structure 

In this chapter, I consider a realistic pillar- like multilayered structure used for current in- 
duced reversal of a magnetic layer f64ll65ll66l consisting of a thick ferromagnetic layer, 
whose primary role is to polarize the current, and whose magnetization is pinned, a thin 
ferromagnetic layer that is to be switched, a nonmagnetic spacer layer so that there is no 
interlay er exchange coupling between the thick and thin layers, and nonmagnetic layer 
or lead in the back of the thin magnetic layer. Fig. 16. 11 The lead in the back of the thick 
magnetic layer does not have to be considered, because the thickness of the left ferro- 
magnetic layer is taken to exceed the spin diffusion length Xsdi in this layer; the spin 
polarization of the current is dictated by the thick magnetic layer. As in the symmetrical 
structure considered in the previous chapter, the thickness of the nonmagnetic spacer 
layer tsp is taken to be small compared to the spin-diffusion length in the nonmagnetic 
metal (Af^^ ~ 600 nm), so that both spin accumulation and spin current are constant in 
the spacer, and all the results are obtained in the approximation that tsp = 0. The major- 
ity of the results for the spin-current, spin-accumulation, spin-torque, and effective field 
acting on the local magnetization in the thin magnetic layer are obtained numerically 
using the Mathematical^ software, W7\ and presented in Sec. 16.11 Next, I derive an 
analytical expression for the transverse spin-current at the interface between two ferro- 
magnetic layers within an approximation that the thickness of the thin FM layer tp is 
much smaller than the spin-flip length A^/, but much larger than the characteristic length 
scale of the transverse spin accumulation Aj. Finally, I study the angular dependence of 
the resistance of the structure described above. 



83 




Figure 6. 1 : Multilayered pillar-like structure used for current induced reversal of a mag- 
netic layer. FM2 is a thick ferromagnetic layer with the thickness exceeding and 

(2) 

local magnetization M)^ = cos Oe^ — sin 6ey, Sp is a thin nonmagnetic spacer, FMl is 
a thin ferromagnetic layer with the thickness tp and local magnetization M^^-* = e^, and 
NM is a nonmagnetic back layer. 



6.1 Numerical results 

In order to find the spin accumulation and spin current distributions in the structure de- 
picted in Fig l6.11 1 solve Eqs. (15.91) and (15.101) in the thick ferromagnetic layer, thin fer- 
romagnetic layer, and nonmagnetic back layer, and find the spin current using Eq. (15.51) 
(see appendix IB. 21) . In Appendices IB . 1 1 and IB .21 1 derive the boundary conditions on the 
accumulation and current at the interface between the thick and thin ferromagnetic lay- 
ers, and at the interface between the thin magnetic layer and the nonmagnetic layer; with 
these one can determine the spin current across the entire structure, and consequently 
the spin torque and effective field acting on the thin magnetic layer. In the absence of 
specular scattering at the interfaces, spin current is continuous across the whole system; 
the spin accumulation is discontinuous across the diffusive interfaces, the discontinuity 
being proportional to the interface resistance, and continuous across the interface with 
zero resistance. I assume that there is neither specular nor diffuse scattering at the in- 
terface between the thin ferromagnetic layer and nonmagnetic layer, and that there is no 
specular scattering at the interface between ferromagnetic layers. I obtain all the results 
both in the absence and in the presence of the diffuse scattering at this interface. With- 
out further simplifications I am unable to give analytic expressions for the accumulation 
and current across the multilayer, and I present the numerical results for these quantities 
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as well as the torque and field they create. In all the plots below, the diffusion constant 
is taken to be 10~'^ m^/s in the magnetic layers, and 5-10~^ m^/s in the nonmagnetic 
layers, and the spin diffusion length in the nonmagnetic layers is taken to be 600 
nm. The amount of diffuse scattering at the interface ARj, the spin polarization param- 
eters for the conductivity 7 at the interface and /5 in the bulk of magnetic layers /3 are 
taken from the experimental data on CPP-MR, ARi{l - 7^) = 0.5 fOm-m^, 7 = 0.75, 
P = 0.5. ll23l Other spin polarization parameters are obtained using the results of the 
ab-initio calculations of the density of states for up and down electrons. 11611 

In Figs. I6.2H6.5I I show the total (not per unit length as in the previous chapter for 
a symmetrical structure) spin torque a sin 9tp and effective field hsinOtp as a function 
of the thickness of the thin magnetic layer tp which is being switched for different 
combinations of Aj, Xgdu ARj, and for different angles 9 between the magnetization 
directions in the layers. Two values of Aj are chosen for illustration, Aj = 4 nm, which 
is comparable to the mean free path Xmfp = 6 nm, and a smaller value of Aj = 1 nm; 
Xsdi is taken to be 60 nm and 30 nm. One can see that while the spin torque rapidly 
increases for small but finite tp ^ \j and then gradually levels off (Figs. I6.2H6.3I) . 
the effective field is largest about tp ^ 0.5Aj and then decreases toward zero with tp 
(Figs. 164116.51) : this can be understood as follows. When the thickness of the thin layer 
tp is much smaller than Aj, the spin accumulation in the thin layer is the same as that of 
the thick layer at the interface, and its direction is parallel to that of the magnetization of 
the thick layer M^^-* (note that I do not consider torques created directly at the interface 
between the ferromagnetic layers); therefore, as follows from the definition of a and b, 
Jm± = aM^^^ X M^^^ + bM^^^ x (mJ^ x M^^^) (Eq. 03), only the effective field 
b can exist, since only the second term has a component along M^^\ As tp increases, 
the spin accumulation in the thin magnetic layer rotates away from M^^-* and develops 
a transverse component, parallel to M^^^ x M^^^ i.e., a spin torque a develops. When 
tp becomes larger than Aj, the spin accumulation is further rotated in the thin layer, 
out of the plane of M^^^'' and M.^J\ and thus the component of the spin accumulation 
in the plane of the magnetizations decreases rapidly, i.e., the effective field diminishes 
faster than the torque. As tp increases further there are no additional contributions to 
either the field and torque because they represent effects that are centered at the interface 
with the spacer layer and averaged over the entire thickness of the thin magnetic layer. 
Although the effective field is negligible compared to the torque in the limit of large 
tp, it is noteworthy that at its maximum the field b is at least as large as a, so that both 
spin torque and effective field should enter the equation of motion of the background 
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magnetization, Eq. (I5.8I ). One also notes that while the torque and field terms a and b 
are largest for 9 = 180° and 0°, they do not act on the background magnetization because 
sin 6' = 0. The largest effects are found for 9 close to 180° (see below). Comparing the 
figures I01 and l6.3ll6.4l and l6.51 one can see that interface resistance decreases both spin 
torque and effective field, and causes the spin torque to achieve saturation for slightly 
smaller tiT. Comparing the diagrams a) and c), b) and d) at each of the Figs. l6.2ll6.51 one 
notices that when the spin diffusion length is reduced from Xsdi = 60 nm to Xsdi = 30 
nm, the spin torque and effective field are reduced. Comparing the diagrams a) and b), 
c) and d) at each of the Figs. l6.2H6.51 one can see that when the spin transfer length is 
reduced from A j = 4 nm to Aj = 1 nm, the spin torque achieves saturation for smaller 
tp and the effective field is increased and peaks for lower t^?. 

In Figs. I6.6|I6.9I I show the spin accumulation and spin current distribution in the 
thick FM-thin FM-NM system both in the absence and in the presence of the interface 
resistance for Aj = 4 nm, Xsdi = 60 nm, for different angles 9 between magnetization 
directions in the layers, and for two different thicknesses tp of the thin ferromagnetic 
layer. I choose tp = 3 nm, smaller than Aj, but where the spin torque a sin 9tp starts to 
saturate, and tp = 15 nm, much larger than A j, but smaller than Xgdi ■ The magnetization 
in the thin layer M^^'' is taken as the global z axis; the current is along the x axis which 
is along the growth direction of the multilayer, and the y direction is perpendicular to the 
other two. In these global axes, longitudinal and transverse in the thin magnetic layer 
refers to the directions z and x — y; however for the thick layer, whose magnetization 
M)^^-* is at an angle 9 relative to M^^-*, the global y and z axes do not define what is meant 
by longitudinal and transverse in this layer. 

Comparing the figures l6!6l and l6.8ll6.7l and l6.91 one notices that the diffuse scattering 
at the interface between the ferromagnetic layers produces sizable discontinuities in the 
accumulation at the interface (see Appendices IB . 1 1 and IB . 2I ). but it does not lead to an 
increased spin torque (see Figs. l6.2H6.5l ): this result is different from the ones obtained 
within other models. 

One can see that far from the interface x <^ —Xj the accumulation and current in 

(2) 

the thick layer are coUinear with background magnetization , i.e., referred to its 
local axes they are longitudinal with no transverse components, and the spin current 
approaches its bare value —>■ /3jeM|^^^ (see Eq. (15. 5I )). as one expects in the bulk of 
a ferromagnet. Even though one still has a longitudinal spin accumulation in the region 
—Xsdi <^ X <ti — A J its gradient is small compared to that of the transverse accumulation 
which makes large contributions to the spin current in the region x > —Xj. This is clear 
from the plots for nix and jm,x which go to zero, while niy — > m sin 9, j^^y (3je sin 9, 
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e = 180° 
e = 150° 
e = 120 ° 

6 = 90° 




9 = 60° 

9 = 30° 
6 = 0° 




Figure 6.2: Total torque asmOtp/ (^jeihal/ 6^3) acting on the thin FM layer in the FM- 
Sp-FM-NM system as a function of the layer thickness for different A j and Xg^i and zero 
interface resistance. 
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e= 180° 

9= 150° 
9= 120° 
9 = 90° 



9 = 60° 
9 = 30° 




Figure 6.3: Total torque a sin OtF/Pjei^CLo/efj.B) acting on the thin FM layer in the FM- 
Sp-FM-NM system as a function of the layer thickness for different Aj and Xsdi and 
non-zero interface resistance. 
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e= 180° 
e = 150° 
e = 120° 

6 = 90° 



6 = 60° 
6 = 30° 
6 = 0° 




tF,nm tp,nm 

Figure 6.4: Total effective field hsinOtp / I3je{ha^/ efJLB) acting on the thin FM layer in 
the FM-Sp-FM-NM system as a function of the layer thickness for different A j and Xsdi 
and zero interface resistance. 
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6 = 150° 

e = 120 ° 

6 = 90° 



O — UL» 

9 = 30° 
9 = 0° 



b sinB tp 
ARj >0 

A. J = 1 nm, X sdl - 30 nm 



a) 



b sinB tp 
ARj >0 

A, j= 1 nm, A. sdl = 60 nm 



c) 




4 6 
b sinG tp 
ARj >0 
Aj=4nm,Asdr 



tF,nm 
d) 



60 nm 




tp,nm 



tp, nm 



Figure 6.5: Total effective field bsinOtp/ Pje{hal/eiXB) acting on the thin FM layer in 
the FM-Sp-FM-NM system as a function of the layer thickness for different A j and Xsdi 
and non-zero interface resistance. 
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TJiz — > mcosO, and j^.z f3jeCos6. For example, for 6 = 90°, jm,y 1, in units 
of Pje, while jm,z ^ in the thick ferromagnetic layer far from the interface. In 
the nonmagnetic layer to the right, x > tp, the spin current is polarized along the 
magnetization direction in the thin magnetic layer 'M.ll\ and there is no, or very little, 
spin current in the x — y direction. 

Surprisingly, there is a huge enhancement of the transverse (x — y) components of 
the spin current in the region about x = (diagrams d) and e) in Figs. l6.6H6.9l) . Since 
the torque transmitted to the thin magnetic layer is the difference between the transverse 
component of the spin current at the boundaries of this layer (see Eq. (IS.MI) '). it is greatly 
amplified compared to what one would find neglecting the spin accumulation. The thick 
magnetic layer to the left x < is pinned so that the enhanced torques acting in the 
region of the interface do not produce any rotation. The z or longitudinal component of 
the incoming spin current is not absorbed by the thin magnetic layer as there is no trans- 
fer of spin angular momentum along this direction {tp <^ ^sdi), see the diagram f) in 
Figs. l6.6ll6.9l The slight decrease in jm,z is due to the spin flip scattering in the magnetic 
layer which is characterized by Xsdi = 60 nm. The much slower decrease in jm,z in the 
nonmagnetic layer, x > tp, comes from the spin flip scattering in the nonmagnetic layer 
whose Xf^i = 600 nm. 

The large enhancement of the transverse spin currents can be understood as follows. 
At the interface between the thick and thin ferromagnetic layers, the spin accumulation 
has to adjust to the new magnetization direction. Both the longitudinal and transverse 
components of m have to experience the change of the same order, but the distance over 
which the transverse component is absorbed, A j, is much smaller than that for the longi- 
tudinal accumulation, Xsdi- Therefore the gradient of the transverse accumulation about 
X = is large and as it is the gradient that contributes to the spin current, Eq. (15. 5I) . one 
finds an amplification of the transverse components of the spin current at the interface. 
In the next section I present a quantitative estimate for the enhancement. 

It is worth noting that since the transport in the spin transfer region is assumed to 
be diffusive (see Chap.[l]), so that the length scale of m± and jm± is Aj, the transverse 
component of the spin current is absorbed over the entire thin layer, if its thickness is 
comparable with Aj (diagrams d) and e) in Figs. l6.6l and l6.8l) . For larger thicknesses tp 
(Figs. l^Jl and l6.9l ). the transverse spin current in the thin magnetic layer indeed goes to 
zero before reaching the interface with the nonmagnetic back layer x = tp.The reason 
for the different behavior in the thin layer arises from its confined geometry, i.e., the 
reflections from the thin magnetic/nonmagnetic back layer interface create the patterns 
observed for the transverse spin currents. 
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Figure 6.6: x-, y-, and 2;-components of the spin-accumulation 
™/(/^Je/"\/2AjJ)(/iao/e;UB) (a-c) and spin-current jm/Pje (d-f) distribution in 
the FM-Sp-FM-NM system with A j = 4 nm, X^^i = 60 nm, the thin FM layer thickness 
tp = 3 nm and zero interface resistance. 
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Figure 6.7: x-, y-, and z-components of the spin-accumulation 
i^/(/5je/v^AjJ)(/iaQ/e/iB) (a-c) and spin-current jm/Pje (d-f) distribution in 
the FM-Sp-FM-NM system with A j = 4 nm, A^^; = 60 nm, the thin FM layer thickness 
tp = 15 nm and zero interface resistance. 
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Figure 6.8: x-, y-, and ^-components of the spin-accumulation 
m/(/5je/v^AjJ)(^aQ/e/iB) (a-c) and spin-current jm/ Pje (d-f) distribution in 
the FM-Sp-FM-NM system with A j = 4 nm, Xg^i = 60 nm, the thin FM layer thickness 
tp = 3 nm and non-zero interface resistance. 
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Figure 6.9: x-, y-, and z-components of the spin-accumulation 
V^^jJ)i^a'o/efiB) (a-c) and spin-current ]m/ Pje (d-f) distribution in 
the FM-Sp-FM-NM system with A j = 4 nm, Xsdi = 60 nm, the thin FM layer thickness 
tp = 15 nm and non-zero interface resistance. 
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6.2 Spin current at the interface between ferromagnetic 
layers 



As seen from Figs. l6.6H6.9l the spin current at the interface between the ferromagnetic 
layers is significantly amplified compared with its bare value ^jj^^ for large angles 9 
between the magnetization directions in the layers. In this section, I derive an analytical 
expression for the amplification within a simplifying approximation \j <^ tp <^ A^/. 

The spin accumulation distribution in the thin layer takes the form (see appendix lB.2l) : 

m^^^ = 2 sin 6e~^ [ReGs cos x — ImG^ sin x] 

+ 2sin^e*-2«[(ReG6Cose-ImG'6sinOcos(x-0 (6.1) 
+ (ImGg cos^ + ReGg sin^) sin(x — ^)], 



m 



(1) 



= 2 sin 9e ^ [ImGs cos x + ReGs sin x] 

+ 2 sin ^e*-2« [(ImGe cos ^ + ReGg sin cos(x - 

— (ReGe cos^ — ImGe sin^) sin(a; — ^)], 



(6.2) 



m 



W = Gse"^/^^* + G4e(^-*^)/^^'^S (6.3) 



where x = x/^/2\j, ^ = tp/^/2\j, G3, G4, G5 and Gq are constants of integration 
found via the boundary conditions (see Appendix lB.il ). The spin current is calculated as 
(see Eq. (1531 ) 

■(1) ^ I gmW(x) 

■Jm,x 2 dx ' 



^.(1) ^ _ldm^{x) 

' ''^ 2 dx 

,■(1) _ ^^n^"^F(a;) 



■(1) ^ y V / .g4x 



2' ' dx 

Using the approximation Xj <^ tp <^ Xsf, the transverse spin-current at the interface 
between ferromagnetic layers x = in the absence of the interface resistance may be 
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31/4 3t/2 331/4 



Figure 6.10: True spin current at the interface between the thin and the thick ferromag- 
netic layers (red line) in comparison with the bare transverse current (blue line) as a 
function of the angle between the magnetization directions in the layers for A j = 4 nm, 
Xsdi = 60 nm. 



written as 



7^^^ (X 



0) = (ReG'5 + ImG'5)sin^/y2Aj, (6.5) 
£,U^ = 0) = -(ReG5-ImG'5)sin^/v^Aj, (6.6) 



or, after substituting G^, 



(3je sin 9 
2(cos2^/2 + Asin^e/2)' 



•(1) / n\ HJe^^f-'^ 



where A = ^1 — /3/3'Aj/a/2As/. Up to the difference in notation, this expression is 
the same as the one obtained in the previous chapter for the spin current at the interface 
between two infinite ferromagnetic layers, Eq. (15.211) . The magnitude of the spin current 
at the interface is shown in the Fig. l6.10l (red line) in comparison with the bare value of 
the transverse current, Pje sin 9 (blue line). One can see that for the angles 9 close to vr, 
the magnitude of the transverse spin current, and, hence, the spin torque acting on the 
thin layer is enhanced compared to its bare value of Pje sin 9, reaching its largest value 

,W_(.^0) = '^^--"';{^ + "> (6.8) 
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for the angle 6* such that 

cosr = -^^. (6.9) 
1 + A 

This huge enhancement, by a factor of 1/4A, where A is of the order of 0.03 for cobalt, 
comes from the interplay between longitudinal and transverse accumulation; it is the 
result of the global nature of the spin current even though the transverse components of 
the spin current and accumulation are absorbed within a region A j of the interface. 



6.3 Resistance 

I calculate the resistance of the structure depicted on Fig. 16. II using the expression for 
the resistivity in each layer which follows from Eq (15.321) for the electric current: 

,6,0, 

2Co jeCo dx 
The sheet resistance of the system AR takes the following form: 

AR = ARi + i^^+ ^-^[Gi + (G4 - G3)(l - e-*-/^-)], (6.11) 

where ARj is the resistance due to the interfaces, L is the total length of the system, Gi, 
G3, and G4 are the constants of integration, defined in the appendix IB. 21 Figure l6?m 
shows the normalized resistance Rnorm = (-R(^) — -R(0))/(-R(7r) — -R(O)) as a func- 
tion of the angle 9 between the magnetization directions in the layers for Aj = 4 nm, 
Kdi = 60 nm, and several values of the thin layer thickness tp. One can see that the nor- 
malized resistance is a non-monotonic function of the angle between the magnetization 
directions in the layers (compare with the reference ll^ ). This is a consequence of the 
angular momentum transfer from the polarized current to the background magnetization, 
and the effect it has on the voltage drop across the multilayer. 
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norm 




3 nm 
15 nm 
60 nm 



Figure 6.11: Normalized resistance Rnorm as a function of the angle 6 between the 
magnetization directions in the layers in the absence of the interface resistance for A j = 
4 nm, Xsdi = 60 nm, and several values of the thin layer thickness tp- 
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Chapter 7 
Conclusion 



Chapters |3l and m address the problem of finding the resistance due to the interfaces in 
the multilayered metallic structures. In Chap.|3l the general method of solving the semi- 
classical linearized Boltzmann equation in CPP geometry is developed, allowing one 
to obtain the integral equations for the chemical potentials everywhere in the multilay- 
ers in the presence of both specular and diffuse scattering at the interfaces, and diffuse 
scattering in the bulk of the layers. In Chap.|4l the numerical results for the chemical po- 
tential and interface resistance in the multilayered metallic systems are presented. The 
variation of the chemical potential within a mean-free path of the interface between two 
semi-infinite metallic layers is investigated. The difference between the potential drop 
associated with the interface /ar from it and that at the interface depends on the amount 
of diffuse scattering at the interface and can change sign. In the case of the same metal 
at both sides of the boundary, the different forms of diffuse scattering are classified 
according to whether the measured resistance is bigger or smaller than that across the 
interface. For certain forms of diffuse scattering, the chemical potential remains con- 
stant outside the immediate interfacial region. Interface resistance in the presence of 
the diffuse scattering in the bulk of the layers is less than that calculated in the ballistic 
transport regime. The resistance due to interfaces in the metallic structures consisting of 
three and five layers is calculated. If one maintains the same interface between layers, 
its contribution to the total resistance of the system depends on the thickness of the adja- 
cent layers. This contribution is constant only if the distance between interfaces is larger 
than the electron mean free path. As the thickness of the layers becomes smaller than 
the electron mean free path in these layers, scattering from interfaces start to interact 
with each other. The error incurred by neglecting the dependence of the interface resis- 
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tance on layer thickness depends on the amount of diffuse scattering at interfaces and 
the height of the potential step at the interface. It reaches 15-20% in Fe-Cr and Co-Cu 
systems with equal amounts of specular and diffuse scattering at interfaces. In order to 
observe this effect one should investigate systems where interface resistance dominates 
total resistance. Hence, the thickness of the layers should be small. Also, by introducing 
additional scatterers, one should be able to change the ratio of thickness to mean free 
path without changing the thickness of the layers. 

In Chaps. |5land|6la mechanism for the magnetization switching that is driven by 
spin-polarized current in noncoUinear magnetic multilayers is studied. I find the spin- 
torque due to the bulk of the magnetic layers and the diffuse scattering at interfaces 
for a number of layered structures; the spin transfer that occurs at interfaces is self- 
consistently determined by embedding it in globally diffusive transport calculations. In 
Chap. 121 a system consisting of two thick ferromagnetic layers divided by a nonmag- 
netic spacer is considered. The spin accumulation and the spin current distributions, 
as well as the spin torques acting on the magnetic layers are found analytically. I ob- 
tain an analytical expression for the CPP magnetoresistance that takes into account the 
effect of the transfer of the spin angular momentum from the polarized current on the 
voltage drop across the multilayer, and estimate a parameter of the angular dependence 
of the CPP- MR that gives a correction to the simple cos^(6'/2) dependence of the re- 
sistance for cobalt and permalloy. For the latter, this parameter can be compared with 
the experimental value. While there is a difference between the theoretical and experi- 
mental values, the theoretical value is obtained without taking into account the interface 
resistance, which presence leads to a smaller discrepancy. In Chap.|6la realistic mag- 
netic multilayer consisting of a thick magnetic layer, a thin magnetic layer that is to be 
switched, a nonmagnetic spacer layer, and nonmagnetic layers or leads on the backs of 
the magnetic layers; the current is perpendicular to the plane of the layers. I find the 
spin-torques acting on the thin magnetic layer, spin accumulation and current profiles in 
the whole structure for different angles between the magnetization directions in the lay- 
ers and different thicknesses of the thin layer. The spin angular momentum transferred 
from the polarized current to the background magnetization of the thin magnetic layer 
far exceeds the transverse component (to the orientation of the magnetization of the thin 
layer) of the bare portion of the incoming spin-polarized current (the part proportional 
to the electric field). Due to the presence of long longitudinal spin-diffusion lengths, the 
longitudinal and transverse components of the spin accumulations become intertwined 
from one layer to the next, leading to a significant amplification of the spin torque with 
respect to the treatments that neglect spin accumulation about the interfaces. I estimate 
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the scale of the amplification as Xsdi/^J, the ratio of the spin-diffusion length to the 
characteristic length scale of the transverse component of the spin accumulation. The 
large enhancment of the spin current in the free ferromagnetic layer may lead to the 
reduction of the critical current necessary to switch spintronics devices. 
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Appendix A 

Appendices for interface resistance 
chapters 



A.l Calculus related to the problem 

In this appendix, several auxiliary equations are obtained, which are used in Chap.|3]to 
derive the equations for the chemical potential profile in multilayers. 



A. 1.1 Dirac delta function in e and and v space 



The Fermi-Dirac equilibrium distribution function has the form 

1 



where e is an electron energy, is Fermi energy, ks is Boltzmann constant, and T is 
temperature. At T ^ 0, the derivative of the distribution function over the energy can 
be approximated by the delta function: 
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The electron's energy is a function of its velocity, so delta functions in e space and v 
space are related via 

e/ ^ N Kv-Vf) 

oe/ov 

and, in a free electron model, 

6{e-eF) = —6{v-VF), (A.l) 
mv 

where m is the electron's mass. 



A. 1.2 Angular averaging 

For an arbitrary function y{ \ cos 6*1), y{ \ cos 6*1) sin6' is symmetric around 9 = n/2, and 
r^%(|cos^|)sin^(i^ = r y{\cose\)smede. (A.2) 

JO Jtt/2 

For example, since a correction to the equilibrium distribution function g^-v, z) can be 
written as a function of z and | cos only, where is an angle that an electron momen- 
tum makes with z axis (Eq. (13.11) ). the chemical potential ijl{z) can be written as 

2/x(z) = rg{e,z)sm9de 
Jo 

= / ((7(1 cos6'|, cos^^ > 0, 2;) + (7(1 cos6'|, cos6' < 0, z)) sin^^d^^ 

JO 



For an arbitrary function y(v), if it can be written as a function of | cos 6'| only, it follows 
from Eq. (lA.ll) and Eq. (IA.2I) that 



/ dv\vz\S(e — eF)y(v) ~ / / v'^dv v\ cos6\-5(v — VF)y(\cos9\) sinOdO 
Jfs Jo Jo v 

= vl r ^ y{\cos9\) cose sin ed9 (A.3) 

JO 
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A.1.3 Diffuse scattering term F{g) in the boundary conditions 

By using the definition of the correction to the equilibrium distribution function g{v, z) 
(Eq. (13.11) ). the term in the boundary conditions responsible for the diffuse scattering at 
an interface F{g) (Eq. (13.71) ) can be written as 

Fi9) = ^[J d^f2\v^M^-eF)g{v<,0) + J d^^,\v>\6{e-eF)g{v>,0)y 
where Vj is an electron velocity in i-th layer, and 

n = J dv2\v<\6{e - ep) + J rfvi|t;> |5(e - e^). 
Then, with the help of Eq. (IA.3I) . 

F{g) = — — / g{vf ,0) cos 9 sin 9d9 

Vpi ~\~ Vp2 Jo 
2i;^ fTT/2 

+ — =^ / c/(f2<,0)cos^sin^ci^ (A.4) 

Vpi + Vlp2 Jo 

is obtained. 



A.2 Current conservation across an interface 

Electrical current j{z) is defined in terms of the distribution function as 

ji{z) = / dvi\vi,\{f{v^,z) - f{v^,z)) 

J FS 

in the left layer, and 

j^[z) = / dvi\v2z\{fiv^,z) - fiv^,z)) 
Jfs 

in the right layer. At an interface, the electron distribution functions f{v{', 0), /(t'f , 0), 
f{v2, 0), and /(f^, 0) are related via boundary conditions (I3.6I) - (I3.7I) . Since the bound- 
ary condition is a linear combination of the specular {S = 1) and diffuse (S = 0) 
scattering terms, it is enough to consider current conservation separately for 5 = 1 and 
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^ = 0. 



If only specular scattering is present at an interface, current to the left of the interface, 
ji{z = 0—, S = 1), may be written as 

= 0-,5 = 1) = 4i r^\oseisme,de,{g{vf,O){l-Ri2{0i))-T2i{e2)9{v<,O)), 
Jo 

and current to the right of the interface, j2{z = 0+), as 

j2(^ = + ,^= 1) =^2 r^%OS^2SinM^2((7K,0)ri2(^l)-(l-i?2l(^2)k«,0)), 

Jo 

where 9i {O2) is the angle of incidence (transmission) to the left (right) of the interface. 
In the process of specular scattering, transmission and reflection coefficients are defined 
so that TijiOi) = 1 — Rij(6i), and angles 9i and 62 are related (Appendix lA. 31) . so that 

Vp^ cos 6*1 sin 9id9i = Vp^ cos 62 sin 6*2 ^6*2, 

and electrical current is conserved at the interface: 

J,(^ = 0-,S=1)=J2(^ = 0+,5=1). 



If only diffuse scattering is present at an interface, current to the left of the interface, 
ji{z = 0—, S = 0), may be written as 

j,{z = 0-,S = 0) = vlJ^'\ose,sme,de^{g{v>,^)-F{g)) 

Jo 

r-K/2 \ 

= vl^ COS 9i sin 9id9ig (vf ,0) - -Vp^F{g) 
Jo 2 

where F{g) is given by Eq. (IA.4I) . Then, expression for ji{z = 0—, 5 = 0) takes the 
form 

Ji{z = 0-,S = 0) = -^^^ cose, sin e,de,{g{v>, 0) - g{v<, 0)). 
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Similarly, current to the right of the interface, j2{z = 0+, S = 0), is 

Uz = 0+,S = 0) = -^^^ r^\ose2sme2de2igiv>,0) - giv<,0)). 

Vpi + Vp2 "'0 

Since angles 6i and 62 are not related in the process of diffuse scattering, it follows that 

3i{z = 0-, 5 = 0) = 32{z = 0+, 5 = 0). 

If both specular and diffuse scattering are present at an interface, so that S* 7^ and 
S ^ 1, the total electrical current ji(2) {z = 0) can be written as a linear combination of 

j^2){z = 0,S = 1) and ji(2)(z = 0,3 = 0): 

Mz = 0-) = Sjiiz = 0-, 5 = 1) + (1 - S)Mz = 0-,S = 0) 

j^{z = 0+) = Sj2{z = 0+, 5 = 1) + (1 - S)j2{z = 0+,S = 0) 

Hence, total electrical current at the interface is conserved by the boundary condi- 
tions (IT6h-(IT7b: 

j,{z = 0-)=j2{z = 0+) (A.5) 



A.3 Specular reflection and transmission at the poten- 
tial step 

Consider an electron with the Fermi energy e p moving in a one-dimensional step-like 
potential V{z) = Vi if z < 0, V{z) = V2 if z > 0. (Tig. lA.lh )) For example, this 
corresponds to a free electron motion across an interface between two layers of different 
metals, or layers of the same magnetic metals with antiparallel magnetization directions. 

If electron energy is referenced from ep, so that = 0, energy conservation law 
takes the following form: 

Vi + = V2 + = ep = 0, (A.6) 

where an effective electron mass m is assumed to be the same in both materials. Electron 
Fermi velocities vpi to the left of the step {z < 0) and vf2 to the right of the step {z > 0) 
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a) 







Vi ~ v2 
^ Fl 






Vo ~ v2 
F2 





£f 
V 



z=0 




z=0 



R2I (62) 



T 12(61) 



Figure A. 1 : Specular scattering at the potential barrier: a) Motion in a step-like potential, 
b) reflection and transmission at an interface. 



are related via 



vfi 

VF2 



\ 



1^2! 



Electron velocity v and momentum k are related as mv = ^k. From the one- 
dimensional geometry of the problem it follows that the component of an electron 
momentum parallel to z = plane is conserved, so that the energy conservation law 
(Eq. (IA.6I )) can be written in terms of z-component of the momentum kz = kp cos 9 
only: 



-2 \- kp^ COS di = -2 V kp2COS 62, 

h n 



(A.7) 



where 9i {62) is the angle that an electron momentum makes with z axis in the left (right) 
layer. 

Condition (IA.7I) means that at an interface, the angle of incidence of an electron 
traveling from left to right 61 and its angle of transmittance 62, or, vice versa, the angle 
of incidence of an electron traveling from right to left 62 and its angle of transmittance 
9i are related as follows (Fig. lA.lb )): 



008 6*2 

sin 62(^62 



cos Oi^ReJl 



sin OidOi^ 



cos^ fl 



+ 



cos-^ 01 



(A.8) 



Vp2 



ReJl- 



1.2 
F2 
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or 



cos6'i 
sin 9id9i 



cos e2^ReJl ^ + 4^^^ 

'^vpi y cos^ 6'2 f cos^ 6'2 

sm^2ci^"^ 



(A.9) 



Re. 1 



cos^ 62 
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Quantum mechanically, an electron is described by the wave function ^p, which, in 
the case of the electron moving in a constant potential, has the form of a plane wave, 
ijj = e*'''". At the potential barrier along z direction, an electron experiences scattering, 
so that its wave function is only partially transmitted across the barrier, and partially 
reflected. For an electron moving from left to right, two plane waves, incident and 
reflected, are present in the left layer, and transmitted plane wave is present in the right 
layer: CFig.lAlb')) 



where r (t) is the reflection (transmission) coefficient for the wave function tp. At an 
interface between layers, both electrical charge density, proportional to and quan- 
tum mechanical flux density, proportional to 'ip'ip*' — ip^'ip', have to be conserved, which 
requires continuity of both tp and ^' at 2; = 0: 

= 0) = = 0), ^[{z = 0) = ij',{z = 0). 

These continuity conditions, together with the expressions for the wave function in the 
left and right layers (lA.lOl) lead to the following expressions for the reflection and trans- 
mission coefficients for the wave function: 




f I ' kiz+k2z 

^2zt I ^ _ 2k 



(A. 11) 



kiz+k2z 



In the boundary conditions (13.61) for the semiclassical Boltzmann equation, reflection 
and transmission coefficients for an electron are required. It follows from the quan- 
tum mechanical flux conservation condition that r and t are related via 

l = |H2 + ^|t|2, 
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so it is convenient to define the reflection (transmission) coefficient for the flux in terms 
of the reflection (transmission) coefficient for the wave function as 



R 



so that the condition T = 1 — i? is satisfied. Reflection coefficient for the electron flux 
going from left to right takes the following form: 



Ri2{ei) = \rM 



1 



1 



I ^F2 1 

+ 7:2 



cos^ 8i v'j^-. cos^ 9i 



1 



1 



"f2 1 



COS^ til Vp, COS^ til 



(A. 12) 



In a similar fashion, reflection coefficient for the electron flux going from the right 
to the left layer can be obtained: 



R 



1 



1 



1 



I ^Fl 1 

+ 7:2 



COS^ 82 v'pr, COS^ 02 



1 -U h 1 I "fi 1 

y ^ COs2 02 ^ COs2 62 



(A. 13) 



Expressions (IA.12I) - (IA.13I) together with the relations on the angles 9i and 62 (IA.9 
A.8I) lead to the microscopic reversibility condition 



R 



R 



(A. 14) 



Figure |A2] shows schematically how the reflection coefficients R12 and R21 depend on 
the cosine of the angle of incidence 61 or 62 for two different geometries of the potential 
barrier, V2 > Vi and V2 <Vi. 



A.4 Equations for chemical potential for three-layered 
system 

In this appendix, equations for the chemical potential profile in the system consisting 
of three metallic layers (Fig. IA.3l) are presented without derivation. The same electron 
relaxation time r is assumed in all three metals, and a dimensionless variable ^ = z/k 
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Ri2 t 

V2 > Vi 1 



R2I 



1 cos 9 I 




R21 



1 cos 9 I 



1 cos 9 2 




1 cos 9 2 



Figure A.2: Reflection coefficient at the potential barrier as a function of the angle of 
incidence 



is introduced, where = vpiT is an electron's mean free path in the i-th layer, so that 
— — ^Ih), and di — d/k, where 2d is the thickness of the middle layer. 



pn/2 

Jo 



(1 + SRi2{ei)) cosOi - ^ST^2{0i) cos ^2(^1) 



X 



exp — sm UidUi 

\ cos 9i ) 



^F2 



(A.15) 
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Figure A. 3: System consisting of three layers. 
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,exp(^^^)sin^i(^2)rf^i(^2) 
cos6'2/ \cos6'i(6'2)/ 



A(^2) exp 



-2^2' 

cos^ 



sin^d^ r^%xp ( ^^^] sin 9d9 
Jo \ COS 9 J 



2/^2(0 = f'\2(e)de exp 

J-di Jo 



+ 
+ 



d2 
7r/2 



COS 6*2 



tan 92d92 



A(^2)expf^-^) sinM^2 

V COS 62 j 



(A. 16) 



/ 5(^2) 
^0 



exp 



-e-^2 
COS 6*2 



sin ^2C?^2 



2/^3(0 



-eE2,h 



rn/2 

Jo 



(1 + SRs2{93)) COS ^3 - T'5r32(^3) COS ^2(^3) 



^F2 
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X exp( ^VM sinM^3 (A. 17) 

V cos6'3 / 



X fexp f-^^") + SR,2{es) exp ( ~^~^'^^'^' ] ] tane.des 

JO \ \ cost's / V ^^^^3 // 

J -d2 

X / ST32{03) exp — — exp — ftftO^^s 

\cos6'2(6'3)y \ cos6'3 J COS62{03) 

X / exp — smd 

Jo \ cos 9 J 

rn/2 



Jo \ cosB J 

+ r^'sT2s{e2)B{e2)exp (^] exp (^^^] sin ^3(^2)c^^3(^2) 
^0 \cos92j \cos93{92) J 

+ — T, 7. — / B\92) exp sm9a9 / exp — sm9a9 

Vp3 + Vp2 JO V ^ / V ^ / 

where 

A{92) = eE2l2\cos92\{l + SR23{92))-ST2s{92)eE3k\cos93{92)\ 

+ SR23{92)B{92) exp (^] + ^^3^ H^' b{9) exp f sin 9d9 

\C0S^2/ ^F3 + ^F2 ■'O \ COS ^ / 

+ si^.3(fe)^exp(^:^)^^^Md^' (a.i8) 

^-da V COS 92 ) COS ^2 

Jds \C0S6'3(6'2)/ C0S6'3(6'2) 
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2 ^ 2 / /^2(e / exp 

Vp^ + •'-d2 Jo \ COS 6 



+ ^'7'' I liz{Odi' r exp ( ' ' ) sin9d9 



Vp^ + Jdi Jo \ COS 9 



B{92) = -eE2l2\cos92\il + SR2ii92)) + ST2ii92)eE,h\cos9,{92) 



SR21 


(^2)^(^2; 






SR21 






J~d2 




r-di 


ST2i{ 


%) / 




J —00 


2(1- 


- SK2 


4i 




2(1- 


- sH, 


4i 


+ VF2 ' 



+ ^ -2^ / ^(^) 1 sin^d^ 



y COS 6*2/ "^Fl + '^F2 "'0 \COs6' 

^^^^^ 

y COS 6^2 / costy2 

^ VcOS^l(^2);COS^l(e2) ^ 
rd2 /-tt/S / —Z^' — r/o \ 

+ ^V^^/ /^2(m'/ exp^-^ sin^^rf^ 
-d2 JO \ cosfc' / 

+ — 2 — : — 2 — / Z^uO"^ / exp —\sm9d9 

-00 JO V cosb' / 



A.5 Equations for chemical potential for five-layered 
system 

In this appendix, equations for the chemical potential profile in the system consisting 
of five metallic layers (Fig. IA.4I) are presented without derivation. For simplicity, only 
symmetrical structures M1-M2-M3-M2-M1 are considered. The same electron relax- 
ation time r is assumed in all three metals, and a dimensionless variable ^ = 2//, is 
introduced, where Zj = Vpir is an electron mean free path in the i-th layer, so that 
f^i = /^i(^ = z/h), and di = d/U, hi = h/U, where 2d is the thickness of the middle 
layer, 26 is the thickness of the second and fourth layers. 

f-di-2bi fTT/2 ( 

2^i(0 = / fiiiOd^' exp tan^ic^^i (A.20) 

J-oo Jo V COSt/i / 
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Figure A.4: System consisting of five layers. 
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2/^3(0 = / l^^{C)di' exp tanM^3 

i-ds "'0 \ COSt^s ) 

+ C{9s) ( exp f - exp f ] sin ^a^^a 



(A.22) 
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where 

A{9i) = e£;iZi(l + SRi2{9i)) cos^i - eE2l2ST^2{0i) cos ^2(^1) 

+ SRuiOl) exp —d^ (A.25) 

\ costal J COS yi 

2(1 - p<il-2''l ^-/2 

+ ,;2 , ,,2 / / exp 

Uz?i fin J —00 jO 



+ -STi2(^i / exp — — — TJK\^^ 

J-d2-2b2 \ cos02{0i) J cos02[0i) 

+ 2 , 2 / i^2{0di' / exp ^ — \ — ^ si^ede 



vpi -r OF2 

+ 5712(^1)52(^2(^1)) exp 



d2-2b2 Jo V COS 6' 

-2b. 



COS ^2(^1), 

+ ^(^^^ r''i?2(^)exp f Z^^l cos^sin^d^ 

■^Fl + '^F2 "'O V COS b* / 



Si(^i) = -e£;2i2(l + 5i?2i(^2))cos^2 + eEi/i5T2i(^2)cos^i(^2) 

+ exp ^< (A.26) 

J-00 \ COsOi[02) J COSUi[U2) 

+ ^ — / /^lU / exp smOde 

Vpi + Vp2 Jo \ COS 6* J 

+ SR2i{e2) exp —d^ 

J-d2-2b2 \ C0S6'2 / C0S6'2 

2(1 - s)vi2 r'^ f -e-d2-2b2 \ . „ 

+ — 2 — ^ — 2 — / A«2(C )"? / exp sm9d9 
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B2{»2) = eB2l2(l + 5^23(92)) cos92-eB3ls5r23(92) 00893(92) 
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A.6 Numerical procedure of solving the Fredholm equa- 
tion of the second kind 

In this Appendix, the numerical procedure of solving the Fredholm equation of the sec- 
ond kind, Eq. (13.41) . will be briefly discussed. The general form of the Fredholm equation 
of the second kind is the following: 

f{x)= f' K{x,y)f{y)dy + g{x), (A.29) 

where f{x) is an unknown function defined on the interval [a, b] of the real axis, K{x, y) 
is the kernel of the integral, and g{x) is the source. 

The Nystrom routine is used in to solve the equation (IA.29I ). It transforms the 
integral equation into the system of N linear equations for the unknown vector f{xi), 
where Xj are the points of the mesh at the interval [a, b], given by the A^-point Gauss- 
Legendre rule. This standard rule ll69l gives the points of the mesh with weights: 

gauleg{a, b, x, w, N), 

where a and b are the boundaries of the interval of interest, x is the vector of the points 
of the mesh, w is the vector of the weights, and N is the number of points in the mesh. 
The equation for the vector f{xi) takes the following form: 

N 

fi^i) = J2 ^(^i' yj)fiyj)'^j + 9ixi), (A.30) 

or 

TV 

fi = J2Ki,f,+gi, (A.31) 



118 



where Kij = WjK{xi, yj), or, in the matrix form, 

(i - k)f = g. (A.32) 

In the case of a singularity in the kernel, for instance if K{x, — > oo when x = y, 
this singularity is removed from the integral by writing 

rb rb 

/ K{x,y)f{y)dy = K{x,y){f{y)-f{x))dy + r{x)f{x), 

J a J a 

where r{x) = K{x, y)dy. Instead of the Eq. (IA.31I) . the following system of equa- 
tions is obtained: 

N _ N _ 

/,(l-r,+ K,,)= K,,f,+g,. (A.33) 

The systems of linear equations (IA.31I) or (IA.33I) are solved using the ludcmp and 
lubksb routines. [69| In order to work with the vectors and matrices, the standard NR.h 
package is used, which contains also the standard numerical routines gauleg, ludcmp, 
and lubksb. 

The numerical procedure for solving the system of the Fredholm equations of the 
second kind, Eq. (13.51) . is identical to the one for solving a single equation, and will not 
be covered here. 



A.7 Parameters entering the equations for the chemical 
potentials for Co-Cu and Fe-Cr systems 

A.7.1 Electrons Fermi velocities ratios for Co-Cu system 

In this appendix, the ratio of the electrons Fermi velocities in the Co-Cu structure 
vfcu/vfco for majority and minority electrons is obtained using the results presented in 
Ref. lITm for the transmission probabilities averaged over the Fermi surface. The proba- 
bility for an electron to be transmitted through an interface between metals i and j as a 
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function of the angle of incidence takes the following form (see Appendix IA.3I) : 



Tij (cos 9) 



4^1 + (^y^|,-l)/cos2 



(A.34) 



l + ^l + (t;Vt;2,-l)/cos2^ 



where Vp^ and Vpj are the electron Fermi velocities in the metals i and j. After aver- 
aging over the Fermi surface, the expression for the transmission probability takes the 
following form: 



where A = Vpj/vp^ — 1. The equation T{A) = T* is solved numerically, with T* 
equal to 0.44 for the minority electrons, and 0.66 for the majority electrons. |70| In 
this illustration, a (110) interface has been chosen, and the electrons were assumed to 
go from Cu into Co. The ratio vpcu/vpco is found to be 1.05 for majority electrons, 
and 1.18 for minority electrons. If the motion from Co into Cu is considered, the same 
ratio of vpcu/vpco is obtained for majority electrons, but different ratio is obtained for 
minority electrons; this suggests that the free electron approximation used to obtain the 
transmission probability expression (IA.34I ) is reasonable for majority states, but not for 
minority states. Within the free electron model, the averaging over the Fermi surface is 
performed, and the information about the true electron momenta is lost, which leads to 
different values of vpcu/vpco for different directions of motion. 



A.7.2 The summary of the Fermi velocities ratios entering the equa- 
tions for the chemical potentials 

The following table summarizes the correspondence between the equations that have 
to be solved in order to find the resistances of up and down electrons R t and -R | in 
different geometries and magnetization configurations (Fig. 14.91) and the electron Fermi 
velocities ratios entering these equations: 




(A.35) 



1 + y^l + 
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Two layers 


Three or five layers 


Eqs.(l3.14ll-(l3.15ll 


Eqs.(IA.15l)-(IA.19l) or (IA.2()D-(IA.28D 




parallel 


anti-parallel 


CoCu 


R T 


^ = 1.05 

Dpi 


fF2 UF2 1 05 


2£z = 1.05; ^ = 1.18 

Vp-l Vpo. 






'"F2 — VF2 = X 18 
Vpx Vp3 


^ = 1.18; ^ = 1.05 

Vpi Vp3 


FeCr 


R T 


^ = 0.837 

DPI 


^F2 — VF2 — gg"^ 

DFl Vp:^ 


^ = 0.837; ^ = 1.003 

VP^ DF.'I 




^ = 1.003 


VF2^Vp2= 1.003 
vp^ vp^ 


^ = 1.003; ^ = 0.837 

VPI Vp^ 
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Appendix B 



Appendices for torques chapters 



B.l Derivation of the boundary conditions for spin-ac- 
cumulation and spin-current between the layers 

In this appendix, the boundary conditions at the interfaces between the layers in Fig. ll.l II 
are derived. To achieve this goal, I consider a sub-system shown in Fig. IB. II which 
consists of a semi-infinite FM layer with a; < with the local magnetization M = 
cos — sin 9ey, a diffuse interfacial layer I between < x < with the same local 
magnetization as in the FM layer, and a semi-infinite NM layer for x > . When 
d^ is infinitesimally small, this sub-system represents the three FM-NM interfaces in 
Fig. 11.1 II i.e., between the thick FM and spacer layers, between the spacer and thin 
FM layers (when spatially inverted), and between the thin FM and back NM layers. 
Both spin accumulation and current are assumed to be continuous at the FM-I and I- 
NM interfaces, and one can derive the relation between spin accumulation and current 
at X = with the same quantities at x = d^ as the thickness of the interfacial layer d^ 
goes to zero. In this limit the parameters of the interfacial layer, such as A^^^, r/^, J^, 
A J, and, most important, its resistance ARj remain constant; the latter condition implies 
that the diffusion constant of the interfacial layer £)g ~ rf^ as rf^ ^ 0. 

Equations ( 15.91 ). (15.101) for the spin accumulation are solved, and the Eq. (15.51) is used 
to find spin current in the ferromagnetic and interfacial layers. By adopting a set of local 
coordinates (x, y, z) such that the local magnetization is Mj = the spin accumulation 
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y 



FM 



I 



NM 



X 







Figure B.l: Structure of the interface between the layers at Fig ll.l II FM is a semi- 

(2) 

infinite ferromagnetic layer with the local magnetization = 0086*62 — sinOey, I is 
a diffuse interfacial layer with the same local magnetization as in FM layer, and NM is 
a semi-infinite nonmagnetic layer. 



and current in the FM layer take the form: 




(B.l) 



and 




(B.2) 
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In the interfacial layer, spin accumulation and spin current take the form: 



= 2i?e(G5exp(^)) + 2i?e ( exp(-f 



' = 2/m G5exp(yf) ) + 2/m G6exp(-y 



ml = Gsexp{^) + G^exp{-^), 

sdl sdl 



(B.3) 



and 



im,y 



Re ( f^exp(^ 



Re ( pi^exp( 



Im ( exp(^) ] - Im[j^ exp(- 



sdl 



G3exp(3^) - G'4exp(-3^ 



sdl 



(B.4) 



Here /?, (3' are spin-polarization parameters for the conductivity and diffusion constant 
in the bulk of the ferromagnetic layer, defined in Chap.|5l 7, 7' are similar parameters 
for the conductivity and diffusion constant in the interfacial layer, defined as 



Ci = 7CoMd, 



(B.5) 



and 



{I 



\ 



l-i 



(B.6) 



when Aj'^ < Af/. 



The boundary conditions for the continuity of the spin accumulation and current at 



124 



the interface between ferromagnetic and interfacial layer x = take the form: 

' 2ReG2 = 2ReG5 + 2ReG6 

< 2ImG2 = 2/mG5 + 2ImGe (B.7) 
Gi = G^ + G4, 

and 

'+ '+ 

^ sdZ sdl 

To relate m^(0) to m.\d^), and j^(0) to jml*^^)' oris can use the assumption that 
as the thickness of the interfacial layer goes to zero, other parameters of the interfacial 
layer, such as J^, and Aj remain constant, but the diffusion constant Dq goes to 
zero with the same rate as d^, so that / = const. Then, for example, for small 
d^ -C Aj the x-component of the spin- accumulation at x = d^ may be written as 

mi{d' 0) ^ 2Re{G, + Gg) + 2Re (^{G, - Gq)^^ . (B.9) 



By comparing this expression with Eqs. (IB.7I) and (IB.8I) . the following relation be- 
tween the x-components of spin accumulation and current at x = and x = d^ is 
obtained: ^ 

mi(d' 0) = (0) - J™,.(0)^, (B.IO) 



and similarly, 

_ 0) = (0) - j^,5(0)^, (B.ll) 



'0 



In a manner similar to Eq. ( IB.9I ) the x-component of spin current at x = d^ may be 
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written as 



)) 



/ Tl 



By comparing this expression with Eqs. (IB.7I) and (IB.8I) . one finds the continuity condi- 
tion for the x-component of spin current: 



and, similarly, 



and 



(B.13) 

(B.14) 
(B.15) 



With these relations, the boundary conditions at the three interfaces in the multilayered 
structure depicted in Fig. ll.lll can now be obtained. 



By using the conditions Eqs. (IB.10I) - (IB.15I) . the boundary conditions at the inter- 
face between the thin (first) ferromagnetic and non-magnetic (N) layers of the structure 
shown in Fig. ll.lll at x = tp may be written immediately, since in the thin FM layer the 
local coordinate system (x, y, z) coincides with the global axes (x, z) \ one finds 



and 







= -rj 






= -rj 






= rje^ 



(1) 

m,x 
(1) 



-77 



-77 



it I 



?(^) it, 



d' 



(B.16) 



(B.17) 
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where r = /2Dq. Similarly, the boundary conditions at the interface between the 
non-magnetic spacer (S) and the thin (first) FM layer at a; = take the form: 

' ^f(0)-m«(0) = rj«,(0) 

< m^(0)-m«(0) = rj«^(0) (B.18) 
_ mf (0) - m«(0) = -rjei^ + rj^Mj^, 







= m«(0)^ 


< Jm,2;(0) 




= _^(i)(0)^ 


V 


-£i(o) 


= m«(0)^, 



Note that spin-current conservation condition at the interfaces, which means that 
there are no torques acting at the interfaces, is due to the infinitely small thickness of the 
interfacial layers ^ 0. To write the boundary conditions at the interface between the 
thick FM and NM spacer layers at x = 0, one have to change from the local coordinate 
system (x,y,z), related to the magnetization direction in the thick FM layer, to the 
global (x, z) system. Any vector a will be transformed according to the following 
rule: 

< Gy = ay cos 6 + sin 9 (B.20) 

ttz = —ay sin 6 + az cos 6. 

By applying this transformation to the conditions, Eqs. (IB.10I) - (IB.15I) . one obtain the 
following boundary conditions at the interface between the thick (second) FM and non- 
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magnetic spacer (S) layers: 

^f(0)-m(2)(0) 



and 



m^(0) -mj2)(o) 



- rjiU^V-^f^ 



+'^jm,z (0) sin 6' cos ^ Y 



77 



-77 



77 
-77' 



mf (0) - mf> (0) = rje cos 9 + rj^/^ (0) sin 6 cos 0^ 



j1.(0) -ji%(0) = -mf (O)cos^^ - mf (O)sin^^ 
j1,,(0) - j;l,(0) = mf (0) cos^^ - mf (0) sin^ 

+mf (0) sin ^ cos 



■sf 



jL. (0) - jf UO) = m(2) (0) sin + m2,(0) sin ^ cos 9 



,(2) 



-mf (0) cos2^4 



2nd' 



(B.21) 



(B.22) 



Note that as the thickness of the interfacial layer goes to zero, for diffuse scatter- 
ing considered here, large discontinuities in the spin- accumulation (Eqs. (IB.16I) . (IB. 181) . 
(IB.21I) ') are produced, proportional to finite r = /2Dq, but small discontinuities in 
the spin-currents (Eqs. (IB.17I) . (IB.19I ). (IB. 221) ) proportional to d^, because Jj does not 
increase and r/^ does not decrease as rf^ 0. Since the torque acting at the interface 
is proportional to the discontinuity of the spin-current at the interface (see Eq. (15.141 )'). 
the finite thickness of the interfacial layer is essential for torque production at the inter- 
face. As the infinitely small interfacial thicknesses d^ ^ are considered, one obtains 
spin-current conservation conditions at each interface: 



(B.23) 
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j^(0)=£H0), (B.24) 

and 

£H0)=£(0). (B.25) 

By eliminating m^{0) and j^(0) from Eqs. (IrTbIi . (IrITT) . (lR24ll . and (lR25ll . one 
obtains the boundary conditions at the interface between thick (second) and thin (first) 
FM layers at x = (of course there is the NM spacer in-between, however its thickness 
tAT is irrelevant for these boundary conditions as long ast^ <^ ^^di)'- 

m«(0)-m(2)(0) = -2rj«,(0) 
m«(0) - m(2)(0) = -rj,^ sin^ - rjO^U^ f-'^^^^^f'^ 

+rjW,(0)sin^cos^^ (B.26) 
m«(0) - m(2)(0) = rj;^(l + cos^) + rjW^(O) sin ^ cos 

■(1) /QN 2-77' sinH 

and 

£n0)=£H0) (B.27) 

Finally, I show that parameter r = /2Dq is proportional to the interface resistance 
ARi found from CPP transport measurements ll23ll . By considering the expression 
(I5.3I) for the electrical current in the interfacial layer, and the assumptions that ~ 
and Xl^i remains constant as the thickness of the interfacial layer d^ ^ 0, one finds 

ARj = d^ /2Cq, or r = ^ = ARjj^. Parameters Cq and Dq may be related via 
Einstein's relation 

Cj = e^Ni{eF)Du (B.28) 

where 

' C = Cj + (T-C 

D = Doi + a- ■'D (B.29) 
^ N = Noi + a- ■'N, 

a is the Pauli matrix. In Eqs. ( IB.29I) . as well as in the equations below, the subscript / 
is omitted for simplisity, as only the parameters of the interfacial layer are considered. 
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Using the definitions of 7 and 7', Eqs. ( IB.5I ). and introducing a spin-polarization 
parameter for the density of states 7", such that 



(B.30) 



equations (IB .291 ) may be written as 



C = Co{I + ia,) 



The product of N and D then takes the form 

ND = iVoDo((l + 77")^ + (7 + i>z), 
and from the Einstein's relation (IB.28I) it follows that 



(B.31) 



7Co = e2iVoDo(7' + 7")- 



(B.32) 



The experimental value of 7 may be found in the Ref. ll23ll . the value of 7" may be found 
from the density of states of up and down electrons, N-^ and A^^^ correspondingly, since 
N may be written as 



N 









and comparing this matrix with the expression for N in (IB.31I ). one obtains 



1' 



(B.33) 



From the equation (IB.32I ) the spin-polarization parameter 7' can be expressed in terms 
of 7 and 7" as 



7 



7-7 
1 — 77' 



(B.34) 
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Figure B.2: Three-layered structure used for current induced reversal of a magnetic 
layer. FM2 is a thick ferromagnetic layer with the thickness exceeding Af^^ and local 
magnetization = cos 6*62 — sin^^ej^, FMl is a thin ferromagnetic layer with the 
thickness tp and local magnetization M^^^ = e^, and NM is a nonmagnetic back layer. 

so that the parameter r takes the form: 

r = ARie^N'^^ep)- (B.35) 

1 — 77' 

where e is the electron charge , Nl{ep) is the density of states at the interface at Fermi 
energy. 



B.2 Solution of the diffusion equation for spin-accumu- 
lation 

In this appendix, I present the expressions for the spin-accumulation and spin-current in 
the multilayered system shown on the Fig. IB. 21 

Equations. (15.91 ). (15.101) for the spin accumulation in each of three layers are solved, 
and the spin currents are found using Eq. (15.51) . In the thick ferromagnetic layer, spin 
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accumulation and current take the form 

' m^2) ^2i?e(G'2exp(^)) 

mf^ = 2/m {G2 exp(^-)) cos 6 — Gi exp(^) sin 6 

777,(2) ^ 2/77J, exp(^-)) sin 6 + Gi exp(i^) cos 6, 
and 



(B.36) 



jS. = -mRe (^exp(^) 
= -/5je sin e - ADoIm exp(^)) cos 9 
+^^^^exp(^)sin. 

j(2) = cos e - ADoIm exp(f )) sin ^ - ^^a^F^Gi exp(y^) cos 

(B.37) 

In the thin ferromagnetic layer 



X — tp 

1+ 



and 



m« = 2Re (G5exp(-^)) + 2i?e (G6exp( 
m« = 27m (G5exp(-^)) + 2/m (G'6exp(i 



mW = G'3exp(-3^) +G'4exp(^), 



& = 4Do [i?e exp(-^)) - i?e exp(3^)) 



(B.38) 



J = AD, [im exp(-^)) - Im exp( 



x—tp 

1+ 



(B.39) 



.•(1) 

Jm.z 



2^0(1-/3/30 



G'3exp(-^)-G4exp(^) 
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;-l _ / I L ~ l-» 

non-magnetic layer, 



where = ~ "^^3^7' ^^'^ ^sdi spin-diffusion length in FM layer. In the 



and 



m^ = Aexp(-^-^), (B.40) 



J 



N 



2^0^ A (_^-iF 



CLAexp(-^-^). (B.41) 



To obtain the 12 unknown constants A^, Ay, A^, Gi, ReG2, ImG2, G3, G4, ReG^, 
ImG^, RcGq, ImGe, one can use the boundary conditions (see Appendix lB.il Eqs. (IB. 161 ). 
(E261), (E231), and (ESTl)): 



f <(t^)-m«M = -rj«.(t^) 



and 



77 



(B.42) 



(B.43) 



^(i)(0)-m(2)(0) = -2rjW,(0) 
m«(0) - (0) = -rj,^ sin^ - rj^^ (0) ^-^^;^ 
+rj«(0) sin^^cos^^ 
m«(0) - (0) = rj,^(l + cos^) + rj«^(0) sin ^ cos 0^ 

' Jm,z\^J I-77' ' 

where the parameter r is proportional to the interface resistance Ai?/ , r = ARje'^No(l— 
7"^)/(l — 77'), e is the electron charge, Nq is the density of states at the interface, 7, 
7', 7" are the spin polarization parameters for the conductivity, diffusion constant, and 
density of states at the interfaces (see Appendix lB.il) . The other six boundary conditions 
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come from the conservation of spin current at the interfaces: 

iZitF)=t\tF), (B.44) 

£n0)=£H0). (B.45) 

In order to find the torque a and effective field b acting on the thin ferromagnetic 
layer, one can use the definition of a and b (Eq. (15.121) ): 

Jm'-l^ = aM^P X M^^^ + bM^j^ x (mJ^ x M^^^) = -asin^e^ - bsmOey. (B.46) 

From the other hand, 

Jm^l'^ = Jm^^^e^ + Jm^ej^. (B.47) 

Comparing these expressions, one can see that the torque a and effective field b per unit 
length can be obtained by averaging the x- and y- components of the spin accumulation 
over the thickness ti? of the thin layer. 
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